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1.0 SUMMARY 


An improved panel method for the solution of three dimensional flow about wing 
and wing-body combinations with leading edge vortex separation is presented. 
The method employs a three-dimensional inviscid flpw model in which the 
configuration, the rolled-up vortex sheets, and the wake are presented by 
quadratic doublet distributions. The strength of the singularity 
distribution, as well as shape and position of the vortex spirals, are 
computed in an iterative fashion starting with an assumed initial sheet 
geometry. The method calculates forces, moments, and detail surface pressure 
distributions. Improvements include the implementation of improved panel 
numerics for the purpose of eliminating the highly non-linear effects of ring 
vortices around doublet panel edges, and the development of a least squares 
procedure for damping vortex sheet geometry update instabilities. 

The documentation is divided up into two parts: 

Volume I Theory Document 

Volume II User’s Guide and Programmer's Document 

Volume I contains a complete description of the method. A variety of cases 
generated by the computer program implementing the method are presented. 

These cases are of two types. The first type consists of numerical studies, 
which verify the underlying mathematical assumptions of the method and, 
moreover, show that the results are strongly invariant with respect to such 
user dependent input as wing panel layout, initial sheet shape, sheet rollup, 
etc. The second type consists of cases run for the purpose of comparing 
computed results with experimental data, and these comparisons verify the 
underlying physical assumptions made by the method. 

Volume II contains instructions for the proper set up and input of a problem 
into the computer code. Program input formats and output are described. A 
description of the computer program and its overlay structure is also 
presented. 


2.0 INTRODUCTION 


A computer program has been developed for the solution of the subsonic, 
three-dimensional flow over wing-body configurations with leading-edge vortex 
separation. The program provides capabilities for calculating forces, 
moments, and detailed surface pressures on thin, sharp-edged wings of an 
arbitrary planform. The wing geometry is arbitrary in the sense that leading 
and trailing edges may be curved or kinked and the wing may have arbitrary 
camber and twist as long as in real flow it produces only a single well 
developed vortex system. The numerical methods employs an inviscid flow model 
in which the wing and the rolled-up vortex sheets are represented by 
continuous quadratic doublet sheet distribution. Furthermore, wing thickness 
may be represented by linear source distributions. The Kutta condition is 
imposed along all wing edges, and a zero force condition is imposed on the 
vortex core. An iterative scheme is applied to find the strengths of the 
doublet distributions as well as shape and position of the free vortex sheets 
spirals satisfying the nonlinear boundary conditions of the flow problem. The 
code includes two iterative solution procedures: (i) Quasi-Newton scheme and 

(ii) Least Squares Method. The least squares procedure for damping 
unstabilities was developed to alleviate convergence problems for certain 
cases using the standard Quasi-Newton iterative scheme. 

The computer program is written in the CDC FORTRAN Extended (FTN4) language 
for the CDC Network Operating System (NOS). The program uses overlay 
structures and fourteen disk files which include the standard system files 
INPUT (TAPES) for card reading and OUTPUT (TAPES) for printing. The program 
has been checked out and run on NASA Langley Research Center's CDC CYBER 
series computers. 

This method was originally developed by the Boeing Company under contracts 
NASl-12185 and NASl-13833. In order to upgrade the capability of the method 
and the code, a coordinated effort was launched involving contracts NASl-15169 
and NASl-15275 from the Langley Research Center and work conducted for the 
Boeing Independent Research and Development Program. For purposes of 
completeness, the independent Boeing work is included in this documentation. 

The documentation is divided into two parts: 

Volume I - Theory Document 

Volume II - User's Guide and Programmer's Document 

The Theory Document (bound separately) contains a detailed description of the 
theoretical method. Also included are computed results which verify the 
underlying mathematical assumptions of the method and test theory comparisons 
which verify the underlying physical assumptions made by the method. 

The remainder of this volume, the User's Guide and Programmer's Document, is 
organized as follows. In section 4 a brief description of the method is given 
for completeness. Section 5 provides instructions for the proper setup of 
analysis case. Network definitions and arrangements are discussed. The input 
formats are described followed by two example cases. Useful hints for 
practical use of these instructions are also included. Section 6 describes 
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the output formats. Discussions and examples are provided. Section 7 
describes the computer programs. This concludes with a description of the 
program structure, the overlay program, the file structure, common block 
definition and a linkage map of the programs and subroutines. 


3.0 NOMENCLATURE 


a 

AR 

b 

c 

Cn 

Cp 

F 

-> 

F 

G 

K 

I 

A 

I 

M 

Mco 

A 

n 

n 

N 


P 

P 

Pi 

P2 

-> 

P 

Q 

Qi 

% 


free and fed sheet geometry parameter 
aspect ratio 
local span 
chord 

normal force coefficient 
pressure coefficient 

equations determining singularity parameters 
force vector 

equations determining vortex geometry parameters 
equations penalizing panel twist 
panel width 

unit vector along vortex core or network junction 

number of grid point rows on a network 

free stream Mach number 

surface unit normal vector 

normal vector at panel center 

number of grid point columns on a network 

circular arc parameter 

pressure 

isentropic pressure 
second-order pressure 
field point 
point on boundary B 
nine canonical panel points 
panel center 
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NOMENCLATURE (CONTINUED) 




s 

(u,v,w) 



x,y,z 

a 


parametric coefficients defining H 
local semi span used in Smith solution 
perturbation velocity vector components 
free stream velocity magnitude 
perturbation mass flux vector 
total mass flux vector 

average surface value of total mass flux vector 
unit vector along x-axis 
Cartesian coordinates 
angle of attack 


Y 

Y 
A 
A 
6 
C 
0 
0 
X 
A 

A 

n 

\) 

P 


V 1 - Moo 2 

delta wing semi apex angle 
ratio of specific heats 

jump in quantity across singularity surface or line 

change in quantity from one iteration to the next 

fraction of Newton step 

surface vorticity vector 

vortex system orientation angles 

all vortex systems geometry parameters 

vortex system scale factor 

all singularity parameters 

doublet strength 

normal vector to panel edge 

fed sheet scale factor 

Newton iteration step size limiter 
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NOMENCLATURE (CONCLUDED) 


source strength 

perturbation potential 

gradient of perturbation potential 

gradient operator 

vector cross product 


4.0 DESCRIPTION OF THE METHOD 


For the sake of completeness, a brief description of the method is included in 
this document. 

4.1 Theoretical Model 


The flow model used in the Leading Edge Vortex (LEV) Program is 
illustrated in Figure 1. Flow about a highly swept wing at angle of 
attack separates at the leading edge and forms a spiral vortex. Studies 
(refs. 1,2) of the principal vortex indicate that its shape and strength 
are relatively independent of Reynolds number. This apparent lack of 
viscosity dependence suggests that the flow may be regarded as potential, 
with the free shear layer represented either as a vortex sheet or, 
equivalently, a doublet distribution supporting a discontinuity in 
tangential velocity. Since the position of the vortex sheet is not known 
a-priori, this results in a problem governed by the linear subsonic flow 
differential equation 


^2 0 


XX + 


zz 


= 0 , 


( 1 ) 


where 0 is the perturbation velocity potential and by non-linear boundary 
conditions . 


The essential elements of the present flow model are the configuration 
surfaces, the trailing wake, the sheet emerging from the wing leading edge 
and the tip (we call this the free sheet), and the rolled-up core or 
spiral region fed by the leading edge and tip vortex sheets (we call this 
the fed sheet). 


The following boundary conditions are imposed on these elements: 


0 The configuration surface must be impermeable. 

(W^-n) = 0 (2) 

where W is the average surface value of the total mass flux vector 
and n is the surface unit normal vector. 


0 


The free sheet and wake cannot support a pressure difference and must 
form a stream surface. 


Ac 


P2 


= 0 


(3) 


where ^ ^P 2 is the jump in the second 'order pressure coefficient, 
see Section 6.2.2 for definition of C . Impermeable condition 


ft -W^^O 


(4) 
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• DIFFERENTIAL EQUATION 




• BOUNDARY CONDITIONS 

•WING, BODY: 
IMPERMEABLE 

• WAKE, FREE SHEET: 

IMPERMEABLE 
ZERO PRESSURE JUMP 

• FED SHEET; 

ZERO TOTAL FORCE 

• KUTTA CONDITION 



FIGURE 1 FLOW MODEL 




0 


The fed sheet is an extension of the free sheet and feeds vorticity 
to the vortex core (modeled as a simple line vortex). The boundary 
condition governing fed sheet size and core orientation is that the 
total force induced on the fed sheet and core by the rest of the 
configuration be parallel to the core. 

t ©AF = 0 (5) 

A 

where i is the unit vector along the vortex core and AF is the 
force. 

The size of the fed sheet is chosen initially by experience or from 
the conical flow results of Smith (ref. 3). 

0 Kutta conditions are imposed along the appropriate leading, side, and 
trailing edges of the wing in the presence of free sheets emanating 
from these edges. The Kutta condition is controlled by the 
appropriate edge matching condition. 

A S = 0 

where f is the surface vorticity vector and H is the unit vector 
along the junction. 

The configuration impermeability condition, the free sheet pressure jump 
condition, and the Kutta edge conditions determine the solution of 
singularity strengths. The free sheet impermeability condition and the 
fed sheet zero force condition will determine the free and fed sheet 
positions . 

In subsonic flow, compressibility is accounted for by use of the Goethert 
rule which is used to transform the problem into the equivalent 
incompressible problem for solution. 

4.2 Numerical Procedure 


This problem can be represent by the proper distribution of logically 
independent paneling networks, which satisfy either Neumann (analysis) or 
Dirichlet (design) boundary conditions. Shown in Figure 2 is a typical 
paneling scheme for a wing-body configuration. Hyperboloidal 
(Hyperbol ic-paraboloid) panels are used to ensure surface continuity. A 
continuous quadratic doublet distribution is used on the midplane to 
represent wing, wake, free and fed sheet networks. A linear surface 
source distribution can be used to represent the body and wing thickness 
if desired. 

The main features of the numerical discretization and computational scheme 
are: 

1) Geometry input for a network consists of a rectangular array of 

corner point coordinates. These corner points are fitted exactly by 
hyperbolic paraboloid patches (hyperboloidal panels). These exact 
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FIGURE 2 PANEL MODEL 


1 





fits ensure surface continuity. 

2) Discrete values of singularity strength are assigned to certain 
standard points on each network. A local distribution on surface 
singularity strength is obtained by fitting a linear source or 
quadratic doublet form to those discrete values in an immediate 
neighborhood by the method of least squares. An analysis type 
network is employed on the wing (geometry of the wing is specified), 
and a design type network of doublets simulates the free sheet 
(unknown free sheet geometry, zero pressure jump specified). In 
order to insure continuity of doublet strength between panels and 
networks, nine degree of freedom splines are used to describe the 
quadratic panel distributions. 

3) Certain standard points on each network are assigned as control 
points, where boundary conditions are specified. These points 
include panel center points as well as edge abutment downwash points 
in the case of doublet networks. The latter serve to impose standard 
aerodynamic edge conditions automatically (e.g., the Kutta condition, 
zero potential jump at thin edges, continuity of singularity strength 
across abutting networks), in order to produce logical independence 
for each network. The number of boundary conditions on each network 
coincides with the number of assigned surface singularity parameters. 

4) The induced potential and velocity integrals of the influence 
coefficient equations are all evaluated in closed form, although 
standard far field expansions are employed when the control point is 
sufficiently distant from the influencing panel. 

Since the problem is non-linear, an iterative procedure must be used for 
solution. An initial guess must be made for the free and fed sheet 
position. Normally results from Smith's conical flow method are used for 
the initial guess, but the user can also input his own geometry. During 
the iterative solution the position and size of the free and fed sheet are 
updated until all the boundary conditions are satisfied. The standard 
free and fed sheet kinematics which allow this updating are shown in 
Figure 3. A cut normal to the longitudinal axis is shown. The wing 
panels, of course, remain fixed. The angle 0 (theta) associated with the 
free sheet segments are free to change with the exception of the angle 
between the horizontal and fed sheets. The length of the free and fed 
sheet segments are controlled by the parameters A (lambda) and the length 
of the fed sheet segments are further controlled by the parameters v 
(nu). These parameters as well as the panel singularity strengths M (mu) 
are all updated simultaneously using a Newton correction scheme. 

The above vortex system kinematics is, of course, only one of many 
possibilities. A good alternative is the kinematics of Smith (ref. 3) 
shown in Figure 4. Here, in contrast to the standard kinematics, angles 
0 are fixed and lengths £ and core location are chosen as the free 
parameters. Both kinematics schemes will lead to the same converged 
solution. Preliminary studies indicate that Smith's kinematics may 
results in faster convergence. 
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FREE PARAMETERS: 0 v 02- ®4- 

FIXED PARAMETERS: «i, 82- ®3. *4- 05- 8 


FIGURE 3 FREE/FED SHEET KINEMATICS 




FREE PARAMETERS: «2, «3, 84, 85. Yc, Zg 
FIXED PARAMETERS: & 2 - ®3' ®4* ®5 


FIGURE if SMITH'S FREE/FED SHEET KINEMATICS 


4.3 Solution Procedure 


The boundary value problem of wings with leading edge vortex separation is 
nonlinear due to the fact that the shape of the free vortex sheet as well 
as its strength are unknown. The solution procedure must therefore be 
iterative. Two solution procedures are available in the LEV code, ITFLOW 
and LSFLOW. 


4.3.1 Quasi -Newton Scheme, ITFLOW 


The standard procedure ITFLOW uses a Quasi-Newton scheme for the iterative 
solution of the flow problem. The incompressible boundary conditions as 
derived from the compressible formulation by application of the Goethert 
rule, can be written symbolically in terms of the following equations 


F(A,0) = 


G (A,©) = 


(W^ •n) = 0 



=0 


(W^- n) = 0 

A , -► 

I ® AF = 0 


Wing-body 

Free sheet and wake 
Kutta condition 
Free sheet 
Fed sheet 


( 6 ) 

(7) 


where A denotes all the singularity parameters and © denotes all 
the geometric degrees of freedom. The function F symbolizes the 
impermeable boundary condition of the wing and body, equation (2), zero 
pressure jump across the free sheet and wake, equation (3), and the Kutta 
condition. The function G represents the stream surface boundary 
condition of the free sheet, equation (4), and the global boundary 
condition of zero net force acting on the fed sheet and the line vortex, 
equation (5). 


Starting with an assumed initial geometry (i.e., a given set of 
parameters © ), the initial singularity strength parameters A are 
obtained using the set of equations (6) in which ACp^ has been replaced 
by the linear form of the pressure equation (see section 6.2.2). 


To obtain a solution, two phases of iterative procedure are performed 
alternatively. The first phase, which is called subiteration, merely 
produces convergence to the nonlinear ACp^ equation associated with the 
pressure jump boundary condition on the free sheet. The spatial location 
of the free sheets is not updated and the aerodynamic influence 
coefficients remain the same throughout the iteration. The Jacobian 
matrix consisting of only the small perturbation of the functions F due to 
the singularity strength parameters (6 f/6a) can be easily calculated. 


6f 

6A 


AA = -pF 


( 8 ) 
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F is known and denotes the error residual in the satisfaction of the 
boundary conditions of equation (6) at intermediate steps in the iteration 
cycle. P represents symbolically the step size scaling parameter 6 
which is a positive number less than 1 and is chosen small enough (by 
the code) to ensure a decrease in F. Newton's method with this controlled 
step size is used and convergence is usually achieved in 2 or 3 iterations. 

For the second phase, the boundary conditions that the free sheet form a 
streamsheet, and the zero force condition on the fed sheet are 
introduced. In general the initial guess, O , will not be correct and a 
full iteration procedure will begin in which the free and fed sheet 
geometry will be updated. This will require the recalculation of those 
aerodynamic influence coefficients affected by the perturbation of the 
free and fed sheet geometry. 

Small perturbations of equations (6) and (7) from the initial "starting 
solution" result in a set of linear equations governing the perturbation 

(;;}-(:) - 

As in equation (8) F and G are known and denote the error residual in 
satisfaction of the boundary conditions at intermediate cycles. These 
equations are solved iteratively by a Quasi-Newton method with controlled 
step size (see Appendix G of Volume I). The calculation of a complete 
Jacobian (left hand side matrix) which includes the effect of the 
perturbation of geometry, ® , is quite expensive. A new Jacobian is 
computed after every three iterations in the iterative process. Five to 
six iterations are generally sufficient to obtain convergence. 

The convergence history of a typical solution is illustrated in Figures 5, 

6 and 7. Figure 5 illustrates the normal force and residual history. The 
subiteration is now shown. Once convergence for the subiteration is 
achieved the complete boundary conditions are introduced and the full 
iteration begins. The solution should not be considered complete until 
the residual is less than 10 . The case shown had a particularly slow 
convergence with the Jacobian update being made only every 5 iterations. 
More typical cases tend to converge in 5 to 6 iterations with Jacobian 
updates occurring every 3 iterations. Figure 6 shows the progress of the 
free sheet geometry at one station during the iteration. Figure 7 shows 
the corresponding pressure distribution. 

4.3.2 Least Squares Method, LSFLOW 

An alternate iteration procedure is also available for those cases for 
which the Quasi-Newton scheme, ITFLOW, fails to converge. In these cases 
local flow anomalies on the free sheet may cause instabilities which 
destroy convergence everywhere in the solution. These instabilities cause 
excessive panel twist which propagates throughout the free sheet. 
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DELTA WING WITH NEAR WAKE 



2 1 0 -1 -2 -3 -4 -5 -6 -7 


LOG (SUM OF SQUARES OF RESIDUALS) 


FIGURE 5 CONVERGENCE CHARACTERISTICS - RESIDUALS 
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A X/C = 1.0 







One of the simplest methods of damping this instability whenever it arises 
is to limit excessive panel twist. This leads to an additional equation 
that all free sheet panels be untwisted (flat). 


K{ 0 ) 


"•^st 


= 0 


( 10 ) 


where n = Q (X) and Q , and Q . are hyperboloidal panel defining 
quantities.^ Equation (IO 7 combined with equations ( 6 ) and (7) creates an 
overdetermined system of equations for A (singularity parameters) 
and 0 (geometric degrees of freedom). 


The system is solved in a least squares sense after suitable normalization 
to account for dimensional differences as well as desired weighting. 
Equation (10) governing panel twist is not weighted heavily since a free 
sheet made up entirely of flat panels may not in general be a good 
approximation to a stream surface. The instabilities produced by a local 
flow anomaly are severe enough that a very small penalty on panel twist 
force relaxation of the boundary condition causing the local anomaly. 

The procedure for solving the overdetermined equation set is iterative as 
before. At the beginning of an iteration, equation ( 6 ) is solved for A 
as a function of the current 0 using Newton's method with controlled 
step size, i.e.. 


6a 


A A = -pF 


( 11 ) 


This is essentially the subiteration which was discussed previousjy in 
Section 4.3.1. Upon obtaining convergence, a new estimate for 0 is 
calculated by solving the equation 

( 6 C. 6 G 

6 A 6 © 6 © 

■6 © 

in a least square sense, where the Jacobian on the left is evaluated at 
the point A =f( 0 ) as determined from ( 11 ) and 9f/90 is 
calculated from 



6f 

6f 

6f 

6a 

60 

60 


(13) 


We assume here that G and K have been normalized appropriately. 

When using the Least Squares Method, a new Jacobian is computed after 
every two iterations. If cycle of step size reduction exceeds 3 (see 
Appendix G of Volume I), then a new Jacobian will also be formed. 


19 



5.0 USER'S INPUT GUIDE 


In this section instructions are given to enable the user to properly set up a 
flow model and prepare the program input data. Since proper formulation of 
the flow model is paramount in obtaining a solution, considerable description 
of the networks and their characteristics is given. Several examples are 
given illustrating the proper network placement for various configuration 
planforms. A complete listing of the Input Formats and two example cases are 
also included. 

5.1 Capabilities and Restrictions 

The Leading Edge Vortex (LEV) program is a versatile tool for calculating 
flows about a class of configurations with leading edge vortex 
separation. The wing geometry may be arbitrary in the sense that leading 
edge and trailing edge may be curved or kinked and the wing may have 
arbitrary camber and twist. The limiting factor on planform shape is that 
only a single primary vortex system be formed. Configurations for which a 
strong well defined vortex system does not exist in real flow will 
probably encounter convergence difficulties during the solution. This 
includes configurations with less than 60 leading edge sweep, 
configurations with discontinuities in the leading edge which will promote 
the formation of more than one vortex system, and solutions at low angles 
of attack where a well defined vortex has not yet formed in real flow. 

Planforms for which successful solutions have been obtained include delta, 
arrow, and diamond wings with pointed or cropped wing tips, and also 
gothic and ogee planforms. Several of these examples may be found in 
Volume I - Engineering Document, Sections 6 and 7. A variety of camber 
and twists have also been successfully analyzed. Several of these 
examples may be found in reference 4. Again the key requirement in any of 
these solutions is that a single well formed vortex exist in the real 
flow. (Note that it may be possible to obtain a solution on a 
configuration with more than one vortex system on each side of the plane 
of symmetry as long as the systems never coalesce. However, this 
capability has not been explored at the time of this writing). 

The program has a symmetry condition option (NSYMM, card 9) which must be 
set in the input. Normally solutions are obtained assuming a plane of 
symmetry. For asymmetric configurations or configurations at yaw the 
symmetry condition must be defeated and both sides of the configurations 
specified.- (For these cases two vortex systems will be specified). The 
network setup for asymmetric cases will be discussed in section 5.4.6, 
results are shown in Figure 23, section 7.1.2, Volume I - Theory Document. 

An often overlooked capability of the LEV program is to analyze attached 
flow models. The setup of such models is identical to that of the 
separated case except that the free and fed sheet networks are deleted. 

Use of this option allows direct comparisons between solutions that assume 
attached or separated flow. An example of this type comparison is shown 
in Figure 28b in Volume I - Theory Document. 
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The program is valid only for subsonic Mach numbers. The Gothert rule is 
applied to transform the problem to the equivalent incompressible case for 
solution. The flow model size restrictions are given in Table 1 which 
appears in section 5.5. Restrictions are given both for the Quasi-Newton 
scheme (section 4.3.1) and for the Least Squares method (section 4.3.2). 
Note that the number of singularity strength parameters does not 
correspond to the number of panels. This is because in the higher order 
panel method used there is not a one-to-one correspondence between 
singularity unknowns and panels. Control point placement on the various 
types of networks which corresponds to the number of singularity unknowns 
will be illustrated in the next section. Also note that for the 
Quasi-Newton scheme the number of singularity parameters (which could be 
used for an attached flow solution) is greater than the combined number of 
singularity parameters, panel orientation angles, and geometry parameters 
which can be used for a separated flow solution. 

5.2 Network Description 


5.2.1 Networ k Nomenclature 

A network is defined as a portion of the boundary surface on which a 
certain distribution of source or doublet strength is specified, together 
with properly posed analysis (Neumann) or design (Dirichlet) boundary 
conditions. The true surface is assumed to have continuous position, 
slope and curvature. Discontinuities in these quantities are therefore 
limited to network edges. The networks are logically independent in that 
each network contributes as many equations as unknowns to the overall 
boundary value problem, hence networks can be added or dropped without 
total reformulation of the problem. 

Every network is specified by giving the coordinates of an array of grid 
points which is basically quadrilateral as illustrated in Figure 8. That 
is, the array consists of M "rows" or grid points which each contain N 
points. N is the number of columns of grid points. A triangular shaped 
network is achieved by allowing one edge of the quadrilateral collapse 
into a single point. This is accomplished by letting a single grid point 
belong to several rows or columns. 

The sense of M and N defines the orientation of a paneling network. Side 
numbering, corner numbering, grid point indexing and outward direction are 
all defind by the sense of M and N. The vector N corresponds to a column 
of grid points directed in the direction of increasing points, while the 
vector M corresponds to a row of grid points in the direction of 
increasing points. The vector NxM is directed out of the surface. The 
outward sense of a network is important when using source type networks. 
The outward side of a source network must always bound the flow. In 
setting up the geometry for a solution it is also important to know the 
proper side numbering nomenclature. Wake and design type networks such as 
those used for the free and fed sheets and the trailing wake demand a 
specific orientation when being attached to the configuration type 
networks (side 1 must attach). Figure 8 illustrates the proper 
nomenclature for a network. Several data preprocessors are included in 
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the LEV code to aid the user in defining the appropriate network input 
data. Use of these preprocessors will be discussed in a later section. 

5.2.2 Network Paneling 

Generally two types of network paneling may be employed, although other 
arrangements are possible. The two basic types illustrated in Figure 9 
are conical paneling and streamwise paneling. Conical paneling is used 
mainly for wing, free and fed sheet networks, while streamwise paneling is 
used mainly for wake networks. Streamwise paneling may also be used on 
wing networks but may require the use of the more expensive least squares 
method, LSFLOW to obtain a solution. Further discussion on the use of 
streamwise paneling on wing networks is given in sections 5.2 and 6.2 in 
Volume I - Theory Document. 

5.2.3 Network Abutments 


A typical problem will consist of several different networks, representing 
different types of singularity and boundary conditions. Control points 
located at the junction of two doublet networks are assigned to match 
singularity strength across the junction. If only one control point 
exists, doublet value is matched. If there are two opposing control 
points the component of vorticity along the junction is also matched. 

Proper edge matching is dependent on correct abutments between networks. 

In order to ensure correct abutments it is absolutely necessary that 
network paneling match identically along adjacent edges. This means that 
adjacent panel grid points across an abutment must be identical to the 
accuracy of the computer. An additional restriction is that network abut 
along complete edges, i.e., their network corner grid points must 
coincide. Examples of acceptable and unacceptable paneling abutments are 
shown in Figure 10. 

Because of the necessity of achieving proper abutments between networks 
before a valid solution may be obtained, a data check procedure (card 15, 
JDATA CHECK) has been incorporated into the program. It is imperative 
that the data check be performed to confirm proper abutment between 
networks before committing a problem to solution. A discussion of the 
abutment data check output will be given in section 6.1.3, 

5.2.4 Network Types and Uses 


The various network types and their uses are illustrated through the 
following example. The paneling scheme of Figure 2 is schematically shown 
in Figure 11. The network type used for each network is summarized in 
Table 2. Several different singularity types and boundary conditions are 
necessary to properly specify the problem. In the present program eight 
network types are available for modeling a given configuration along with 
its separated vortex system. Each network type represents a different 
source or doublet distribution accompanied by a properly posed set of 
boundary conditions. These network types are distinguished by the index 
NT. A brief description of each available type is presented below. 
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FIGURE 10 ABUTMENT EXAMPLES 
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TYPE AND USE OF EACH NETWORK OF ARROW WING-BODY 

TABLE 2 


Network Sequence Network Use 

Number Type 


1 

NT = 1 

2 

NT = 1 

3 

NT = 1 

4 

NT = 1 

5 

NT = 1 

r 

0 

NT = 1 

7 

NT = 2 

8 

NT = 2 

9 

NT = 8 

10 

NT = 8 

11 

NT = 10 

12 

NT = 6 

13 

NT = 8 

14 

NT = 4 

15 

NT = 14 

16 

NT = 4 

17 

NT = 14 

18 

NT = 16 

19 

NT = 10 


Upper Forebody 
Lower Forebody 
Upper Midbody 
Lower Midbody 
Upper Aftbody 
Lower Aftbody 
Wing 
Wing 

Carry Over Lifting System 
Carry Over Lifting System 
Wake of Carry Over System 
Near Wake 

Trailing Wake of Near Wake 

Free Sheet 

Fed Sheet 

Free Sheet 

Fed Sheet 

Wake of Free Sheet 
Wake of Fed Sheet 
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NT = 1: Source/Analysis Network 

This network is used primarily to represent the exterior surfaces of thick 
wings and bodies. See network numbers 1 to 6 of Figure 11 as examples. 

When inputting source type networks one should always be careful that the 
surface normal (N x M) points out into the flow. 

The singularity parameters and control point locations for the 
Source/Analysis network are illustrated in Figure 12. 
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FIGURE 12 PARAMETER AND CONTROL POINT LOCATION ■ 
SOURCE/ANALYSIS NETWORK 





NT = 2: Doublet /Analysis Network 


This network is used primarily to represent a thin wing and is placed on 
the camber surface of the wing (e.g., networks 7 and 8 of Figure 11). 

This network type is also used as a lifting system for a thick wing. Here 
the network is placed on the camber surface in the same fashion as for a 
thin wing. However, Source/Analysis (NT = 1) networks are then added to 
form the upper and lower wing surfaces. 

The singularity parameter and control point locations for the 
Doublet/Analysis networks are illustrated in Figure 13. 
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FIGURE 13 PARAMETER AND CONTROL POINT LOCATION 
DOUBLET/ANALYSIS NETWORK (NT = 2) 
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NT = 4: Doublet/Design #1 Network 

This network is used as a free sheet, that is, a sheet which has A Cp = 0 
boundary conditions and is updated to be a stream surface. (Cp here is 
calculated using the second order formula, equation 18, section 6.2.2.) 

See networks 14 and 16 of Figure 11 as examples. These examples 
illustrate two important rules concerning the corner point input of a free 
sheet network. First, the apex or collapsed side of a free sheet must be 
side 4. Secondly, the side adjoining the wing (or adjoining another free 
sheet attached to the wing) must be side 1. 

The singularity parameter and control point locations for the 
Doublet/Design #1 network are illustrated in Figure 14. 
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FIGURE 14 PARAMETER AND CONTROL POINT LOCATION 
DOUBLET/DESIGN #1 NETWORK (NT = 4) 





NT = 6: Doublet/Design #2 Network 


M 


FIGURE 


-- 


This network is used for a wake in place of a type NT = 8 network when the 
approximation of the linearized pressure formulas is deemed insufficient 
(see discussion in section 7.1.1, Vol . 1). The boundary condition /^Cp = 
0 (where Cp is calculated using the second order formula, equation 18, 
section 6.2.2), is applied on each panel. In contrast to the type NT = 4 
network, this network must remain fixed. See network 12 of Figure 11 as 
an examaple. If an additional wake is attached to a type 6 network as in 
Figure 11, the wake should adjoin side 3 of the type 6 network. 

The singularity parameters and control point locations for the 
Doublet /Design #2 network are illustrated in Figure 15. 
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NT = 8: Doublet/Wake #1 Network 

This network is used as a wake behind a wing. It satisfies a built-in 
boundary condition, namely that ACp = 0, (where Cp is calculated using 
the linearized pressure formula equation 16, section 6.2.2). This is 
achieved by making doublet strength constant along columns (which are 
presumed to be in the stream direction). See network 13 of Figure 11 as 
examples. One rule concerning corner point inputs for an NT = 8 network 
is that side 1 must always be placed next to the wing or near wake 
trailing edge. 

A type 8 network is also used as a carry-over lifting system which extends 
the wing lifting system into the body (see networks 9 and 10 in Figure 
11). For this purpose the type 8 network is turned sideways. Note that 
side 1 must adjoin the wing lifting system root edge. 

The singularity parameters and control point locations for the 
Doublet/Wake # 1 networks are illustrated in Figure 16. 
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FIGURE 16 PARAMETER AND CONTROL POINT LOCATION 
DOUBLET/WAKE #1 NETWORK (NT ^ 8) 


32 




NT = 10: Doublet/Wake #4 Network 

This network is used as a wake behind a carry-over lifting system or a fed 
sheet. It has constant doublet strength and therefore, carries no shed 
vorticity. As examples see networks 11 and 19 of Figure 11. Note that 
side 1 is always placed next to the carry-over lifting system or fed 
sheet trailing edge. 

The singularity parameter and control point locations for the Doublet/Wake 
#4 network are illustrated in Figure 17. 
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FIGURE 17 PARAMETER AND CONTROL POINT LOCATION 

DOUBLET/WAKE # 4 NETWORK (NT = 10) 
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NT = 14: Fed Sheet Network (Doublet/Wake # 2) 


This network is of the same basic construction as the type NT = 8 network 
but has special panel center and terminated edge velocity evaluation 
points for the calculation of the total force on the network. See 
networks 15 amd 17 of Figure 11 as an example. Note that side 1 must 
adjoin the free sheet. 

The singularity parameters and control point locations for the Fed Sheet 
network are illustrated in Figure 18. 


PARAMETER LOCATION CONTROL POINT LOCATION 




N 


I ^ 

FIGURE 18 PARAMETER AND CONTROL POINT LOCATION 
DOUBLET/WAKE #2 NETWORK (NT = 14) 





NT = 16: Doublet/Wake #.3 Network 


This network is used as a wake behind a free sheet. It is just like a 
Doublet/Wake #1 network (type 8) except that its degrees of freedom are 
associated with its edge corner points so that it can be used behind a 
free sheet. It satisfies a built-in boundary condition, namely that A Cp 
= 0 (where Cp is calculated using the linearized pressure formula, 
equation 16). This is achieved by making doublet strength constant along 
columns (which are presumed to be in the stream direction). See network 
#18 of Figure 11 as an example. Note that side 1 is always placed next 
to the free sheet trailing edge. 

The singularity parameter and control point locations for the Doublet/Wake 
#3 networks are illustrated in Figure 19. 
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FIGURE 19 PARAMETER AND CONTROL POINT LOCATION 
DOUBLET/WAKE #3 NETWORK (NT = 16) 
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5.2.5 Network Geometry Preprocessors 


In order to facilitate the definition of the various networks in the 
program input several geometry data preprocessors have been included in 
the code. These preprocessors will greatly simplify the user's task by 
generating the network grid point geometry for most cases. A preprocessor 
option (^POINTS) is also included to allow completely general grid point 
definition by some outside source. The following preprocessors are 
available within the code. 

^POINTS 

3QUADR I LATERAL 

?G0THIC 

^TRAILING WAKE 

3V0RTEX 

In addition, a sixth preprocesor called ^CAMBERED WING is also available 
for use with 3QUADR I LATERAL or ^GOTHIC for the simplified input of 
cambered or twisted surfaces. 

A description and discussion of the various preprocessors follows. 

^POINTS - This option allows the input of a user defined array of XYZ 
points in 6E10.0 format. It can be used to define any type of network and 
is the most general input for planform and camber shapes. The input 
formats for ^POINTS are given in section 5.5. Users may find some 
difficulty in obtaining proper abutments between networks defined using 
^POINTS and networks defined using ^GOTHIC or ^QUADRILATERAL. This may 
occur because of the requirement for adjacent grid points across abutments 
being identical not being met. The abutments may appear identical in the 
program printout but it must be realized that the grid point coordinates 
calculated by ^GOTHIC or ^QUADRILATERAL may have more significant figures 
than shown in the output format. In such cases it may be convenient to 
dispose the network geometry (which is saved on TAPE14) as punch card 
output, and reinput the networks using ^POINTS. 

^QUADRILATERAL - This option allows the definition of a network by 
specifying the network corner points and the internal percentage arrays to 
define the paneling distribution. This option is useful in defining 
simple wing planforms and design wakes. The input formats for 
^QUADRILATERAL are given in section 5.5. 3CAMBERED WING may be used with 
this option to define a camber and a twist for the network. 

^QUADRILATERAL may be used to generate any type network. 



^GOTHIC - This option allows the definition of a network with a straight 
or curved edge. A longitudinal array of XYZ points is input which defines 
both the edge geometry and the longitudinal panel spacing. A percentage 
array defines a lateral panel spacing. The input formats for ^GOTHIC are 
given in section 5.5. ^CAMBERED WING may be used with this option to 
define a camber and twist for the network. ^GOTHIC may be used to 
generate any type network. 

^TRAILING WAKE - This option is used to define a simple network paneling 
which consists of a single row of panels. This network attaches to 
another specified network edge and extends straight back (parallel to the 
X-axis) to a specified distance. This option is used exclusively to 
generate wake networks. The network and edge to which the wake attaches 
is defined to ensure a proper abutment. The input formats for ^TRAILING 
WAKE are in section 5.5 . 

3V0RTEX - This option will automatically generate a free and fed sheet 
network and their associated trailing wakes. The shapes of the free and 
fed sheet networks are based on Smith's conical results which are 
discussed in section 5.3. The network edge to which the free sheet 
attaches is defined to assure proper edge abutment. More than one set of 
free and fed sheet networks may be specified in tandem in order to satisfy 
the edge abutment constraints when more than one network has been used to 
define the wing. This option may only be used to define free (NT=4) and 
fed (NT=14) sheet (and associated wake) networks. The input formats for 
JVORTEX are given in section 5.5. 

3CAMBERED WING - This option is used in conjunction with 2G0THIC or 
$QUADR I LATERAL to generate network goemetry for cambered or twisted 
surfaces. Camber lines can be defined independent of the network arrays 
to define a 3-D cambered surface. Linear spanwise interpolation is used 
to generate the cambered surface at the network grid points. 

^CAMBERED WING can also be used to generate conical camber of the form 
used by Wentz (ref. 5) or conical camber where the wing is a portion of a 
circular arc, 

5.3 Starting Sol u tion 

The iterative process used in the solution of the leading edge vortex 
problem required an initial guess for the free and fed sheet goemetry. A 
reasonable guess may be based on the conical solution of Smith (ref. 3). 
Smith's results are reproduced in Figure 20, which shows the shape of the 
free sheet and the size of the fed sheet for various values of the 
parameter "a". This parameter is defined as 

® (14) 


where a denotes the angle of attack in radians and *Y is one-half of the 
apex angle of a delta wing. The sheet geometries of Figure 20 represent 
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FIGURE 20 INITIAL FREE-SHEET GEOMETRY AND SIZE OF 
FED SHEET FOR VARIOUS 'a ' 






transverse cuts through the configuration normal to the wing surface. The 
y,z-coordinates are nondimens ional i zed by the wing semi span s. The 
locations of the line vortex along the terminated edge of the fed sheet 
are given for several values of a (0.2 < a ^ 3.0) and are connected by a 
dashline. The straight line between the last point on the free sheet and 
the line vortex is the trace of the fed sheet. An example is shown for 
"a" = 1.4. 

Smith's results are available within the LEV program through a network 
preprocessor called 3V0RTEX. Here Smith's solution in tabulated form is 
used to form the network goemetry for the initial guess on the free and 
fed sheet shape. Figure 21 illustrates how an initial free-sheet geometry 
is obtained for a nonconical wing geometry. For this purpose the 
assumption is made that initially the shape of the free sheet at a 
particular chordwise station is the same as that of a certain delta wing. 
This delta wing is locally equivalent to the considered nonconical wing 
geometry and is defined as a wing that has the same apex position and the 
same local semispan at that chordwise station where the initial free-sheet 
geometry is to be computed. Thus, the parameter a can be calculated at 
each transverse cut for a given angle of attack and a given angle Y - arc 
tan (s/x). Linear interpolation of Smith's data provides the desired 
initial free-sheet geometry for a chosen number of free-sheet panels. All 
free-sheet segments of a transverse cut (y,z-plane) have approximately the 
same chord length. 

The described procedure also provides the size of the fed sheet at all 
geometry defining transverse cuts. During the iteration process the shape 
of the free and fed sheet may change dramatically as shown in Figure 22. 
Here a transverse cut through the free and fed sheet at the trailing edge 
of an aspect ratio 2.0 delta wing is shown. The initial guess is 
generated by the ^VORTEX preprocessor based on Smith's conical solution. 

The converged position shows considerable growth of the free and fed sheet. 

It should be emphasized that Smith's conical data provide only the initial 
free-sheet geometry. This is a convenient choice and a good guess for 
wing geometries that are not too different from flat delta wings. The 
computed doublet distributions and the sheet geometries computed in 
subsequent cycles of the iteration procedure are, in general, not conical. 

The choice of initial and fed sheet shape is in general not critical to 
the converged solution. However, the choice will affect the number of 
iterations necessary for a solution to achieve convergence and may 
preclude convergence in some cases. Figure 23 illustrates a case in which 
an asj^metric initial guess was used. Eight iterations later, the 
solution had converged to a symmetric solution. 

Initial guesses other than those based on the Smith’s results can be 
used. The most general option available is the use of the ^POINTS 
preprocessor. ^POINTS allow a network definition by xyz panel corner 
points, which gives the user complete freedom in specifying free and fed 
sheet shape. Sometimes it is desirable only to grow the Smith guess in 
order to increase the clearance between the fed sheet termination and the 
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FIGURE 23 ASYMMETRIC INITIAL SHEET SHAPE 



wing surface (i.e., on a highly cambered wing). This may be accomplished 
within the 2V0RTEX preprocessor by use of the ARC parameter on card V5. 

ARC is an increment to the "a" parameter of equation 14. A positive value 
for ARC will result in a larger free and fed sheet. In general it is 
easier for a large free sheet /fed sheet to contract to a converged 
position than for a small initial guess to grow into a converged position. 

5.4 Example Network Arrangements 

A series of examples illustrating the proper network arrangement for 
various configurations are presented. The purpose of these examples is to 
aid the user in the proper use and placement of the various types of 
network in formulating a flow model. These examples do not necessarily 
represent the only possible modeling for the various configurations. In 
explaining these various network arrangements it is of particular 
importance to note the network orientation. For those network for which 
the orientation is important, network side number one has been identified 
along with its M and N vector orientation. 

5.4.1 Delta Wing Wi t hout Near Wak e 

A delta wing without a near wake is about the simplest model for which a 
leading edge solution can be obtained. This model is in general not a 
good model because of the inadequacy of the doublet/wake 1 (NT=8) network 
in satisfying the Kutta condition (see section 7.1, Vol. 1 - Theory 
Document). Its use, if at all, should be limited to delta wings of aspect 
ratio less than 1.0 and angles of attack greater than 15 . 

The network formulation schematically is shown in Figure 24 with the 
vortex system rolled out flat in the plane of the wing. The flow model is 
made up of: 1, a doublet /analysis (NT=2) network for the wing; 2, a 

doublet/wake # 1 (NT=8) network for the wake from the wing; 3, a 
doublet/design #1 (NT=4) network for the free sheet; 4, a doublet /design 
#3 (NT=14) network for the fed sheet; 5, a doublet/wake #4 (NT=16) network 
for the wake from the free sheet; and 6, a doublet/wake # 2 (NT=10) network 
for the network for the wake from the fed sheet. 

In setting up this model 2R0INTS, ^QUADRILATERAL, or ^GOTHIC could be used 
to define the wing network, ^TRAILING WAKE to define the wake from the 
wing network and JVORTEX to define the remaining four networks 
representing the free and fed sheets and their associated wakes. 

5.4.2 Delta Wing With Near Wake 

This is the recommended model for most wing planforms. The network 
arrangement is shown in Figure 25. The near wake is actually a 
doublet/design # 2 (NT=6) network, which satisfies the boundary condition 
ACp = 0 where Cp is calculated using the second order formula, equation 
17, section 6.2.2. Unlike the simpler doublet/wake #1 (NT=8) network 
which only satisfies the linear ACp = 0 boundary condition, the design 
wake network accommodates a spanwise shedding of vorticity at the trailing 
edge which is necessary to properly satisfy the Kutta condition. 
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FIGURE 25 DELTA WING WITH NEAR WAKE NETWORK ARRANGEMENT 
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Studies (Section 7.1.1, Vol . I - Theory Document) have shown that the 
design wake can be as short as 0.1 root chords and only two rows of panels 
deep (in the x-direction) . The planform may correspond to a simple 
extension of the wing planform or may increase the sweep of the edge from 
which the free sheet abuts. It is not recommended that the free sheet 
edge become parallel to the x-axis as this can make convergence more 
difficult. The free and fed sheets must extend to the end of the near 
wake. Because of the requirement that abutments occur along complete 
network edges it is necessary to split the free and fed sheets into two 
segments as shown in Figure 25. 

^POINTS, ^QUADRILATERAL, or JGOTHIC may be used to generate the wing and 
near wake geometry. ^TRAILING WAKE will take care of the wake network 
number 3 in the example. $V0RTEX can be used to generate the remaining 
networks, numbers 4, 5, 6, 7, 8 and 9. Two calls to 3V0RTEX will be 
necessary to generate the six networks. First, networks 4 and 5 will be 
generated without trailing wakes (KW=0, Card V2). Then networks 6, 7, 8 
and 9 can be generated. NATF (Card V4) will assure connection of networks 
6 and 7 to networks 4 and 5. JNAT (Card V5) will be referenced to the 
vortex apex network (Network 4) to assure a proper starting solution. 

5.4.3 Arrow Wing- 

A network arrangement for an arrow wing planform is shown in Figure 26. 

The need to split the wing, free and fed sheets into two networks depends 
on the type of paneling used to define the wing. If a conical type 
paneling (see Figure 9, example 2) is used, then the split in the wing, 
free, and fed sheets is necessary as shown in Figure 26 and is due to the 
constraint that networks abut along entire edges (section 5.2.3). If a 
steamwise paneling scheme (see Figure 9, example 1) is used then the wing 
and free and fed sheets could be single networks. 

5.4.4 Rectangular Wing 

A network arrangement for a rectangular wing is illustrated on Figure 27. 
Limited studies (section 7.2 - Vol. I - Theory Document) have not 
indicated a need for a design wake on rectangular wings. Some problems 
have been encountered with the convergence on rectangular wings associated 
with the starting solution generated by JVORTEX when APC=0. (card V5). 
Setting APC = 0.5 to 1.0 helped to avoid the convergence problems. 

5.4.5 Wing With Cropped Tip 

Figure 28 shows a network arrangement for a wing with a cropped wing tip. 
The presence of the tip does not by itself introduce any new network 
arrangement procedures. The arrangement shown assumes that the tip has 
been defined as part of the leading edge so that the leading edge and the 
tip together form only one edge. Conical type paneling would then be 
necessary and could be defined by use of 2G0THIC. Planforms with tips 
cropped parallel to the x-axis have experienced convergence difficulty 
using the standard ITFLOW iteration procedure, necessitating use of the 
more expensive LSFLOW procedure. 
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A simple gambit which tends to alleviate the convergence problems is to 
modify the tip as shown in Figure 28. Instead of a tip with an edge 
parallel to the x-axis, the tip is modified to have a high sweep but still 
retain the same planform area. 

5.4.6 Asymmetric or Yaw Configurations 

The network arrangements for the preceding configuration all assumed a 
plane of symmetry and therefore zero yaw. The paneling arrangement for an 
asymmetric configuration or configuration at yaw demands that the entire 
(both sides) configuration be represented. Such an arrangement is 
illustrated in Figure 29, NSYMM (Card 9) must be set equal to zero. The 
wing may be represented by one of two networks split along the x-axis. 

When two networks are used for the wing, care must be taken to keep the 
upward sense (NxM) of the networks the same. 

5.4.7 Wing Body Configuration 

A network arrangement for a wing-body configuration has already been shown 
in Figure 11. This figure was used as an example to discuss the various 
network types and their uses in section 5.2.3. Table 2 summarizes the use 
of the various networks. One alternative possible over what is shown in 
Figure 11 is to combine all the source networks (1-6) into one network. 
Source networks are exempt from the edge matching requirements of the 
doublet networks. 

5.5 Input Format Specifications 

The input data sequence is illustrated in Figure 30. The data sequences 
consists of several cards defining the flow conditions, configuration 
reference values, program execution mode, etc. These are followed by a 
series of network data blocks which define the flow model. A network data 
block consists of any one of several data preprocessors such as ^POINTS, 
3QUADR I LATERAL, $G0THIC, ^TRAILING WAKE, and ^VORTEX. Two of these 
preprocessors can also include a call to ^CAMBERED WING for simplified 
input of camber surfaces. 

Program size limitations have been summarized in Table 1. A summary of 
the various types of networks and their uses is given in Table 3. 

All numerical inputs are read in 6E10.0 floating point format. Some input 
variables are named in traditional integer format. These designations are 
internal designations and the data should be input as a floating point 
number. All literal words are read in A4 format. Only the first four 
characters need be input. A description of the data input follows. 


50 




FIGURE 29 ASYMMETRIC OR YAW CONFIGURATION NETWORK ARRANGEMENT 









PROGRAM SIZE RESTRICTIONS 
TABLE 1 

Restriction on total number of networks 


NNETT < 20 

Restrictions for Quasi-Newton Scheme, ITFLOW 


NF (No. of singularity strength parameters) < 750 
NG (No. of free sheet panels) < 300 
= No. of panel orientation angles 
NH (No. of fed sheet panels) < 50 
= No. of geometry parameters 
NF + NG + NH < 500 


Restrictions for Least Squares Method, LSFLOW 


NF < 400 

NG + NH < 80 

NG + NH + NK < 144 

Where NK (=NG) is the number of twist function equations 
NF + NG + NH < 480 
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NETWORK TYPES AND THEIR USES 
TABLE 3 



Type 

Common Use 

NT = 1 

Source/Analysis 

Exterior surface of thick wings and 
bodies 

NT = 2 

Doublet/Analysis 

Camber surface of wing 

NT = 4 

Doublet/Design #1 

Free sheet 

NT = 6 

Doublet/Design #2 

Near Wake 

NT = 8 

Doublet/Wake #1 

Simple wake or carry over lifting system 

o 

i-H 

II 

1 — 

Doublet/Wake #4 

Wake behind carry over lifting system 
or fed sheet 

NT = 14 

Doublet/Wake # 2 

Fed sheet 

NT = 16 

Doublet/Wake #3 

Wake behind free sheet 
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Data Card Variable (Numeric data are input in 6E10.0 format) 


Number 

Column Name 

Descr i pt i on /Comment 

1 


3CASE 

2 

1-80 

Title information 

3 

1-80 

User information 

4 


?^NGLE OF ATTACK 

5 

1-10 ALPHA 

Angle of attack in degrees 

6 


gYAW ANGLE 


7 

1-10 

BETA 

Yaw angle in degrees; if NSYMM = 0. both sides of 
the configuration must be defined 

8 



^SYMMETRY 

9 

1-10 

NSYMM 

1. for symmetry about X-Z plane; 0. otherwise 

10 



3MACH NUMBER 

11 

1-10 

AMACH 

Mach number; must be less than 1.0 

12 



^REFERENCES 

13 

1-10 

XREF, 

XREF, YREF, ZREF are the x,y,z coordinates of the 


11-20 

YREF, 

moment center. 


21-30 

ZREF 


14 

1-10 

SREF, 

SREF is the configuration area. 

(Half area if only half is panelled, NSYMM=1.) 


11-20 

BREF, 

BREF = span reference length 


21-30 

CREF, 

CREF = chord reference length 


31-40 

DREF 

DREF = height reference length 


15 


^ITERATION or 2LEAST SQUARES ITERATION or 
JOTITERATION or gPATA CHECK 


IF ^ITERATION or JLEAST SQUARES ITERATION is specified then follow 
with ITMX on next card 


16 1-10 ITMX Maximum number of iterations allowed for the 

iterative procedure. IF ITMX < 0, the program 
will read corner points and singularity strength 
parameters data (from previous run) on disk file 
TAPE14 provided by the user. 

IF ^ITERATION is specified then the program will use Quasi-Newton 
scheme to find an iterative solution of the flow problem. A new 
Jacobian will be computed after every 3 iterations. 
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Data Card Variable (Numeric data are input in 6E10.0 format) 


Number Column Name Descript ion /Comment 


If 3LEAST SQUARES ITERATION is specified then the progrm will use 
Least Squares method to find an iterative solution of the flow 
problem. A new Jacobian is computed after every 2 iterations. 
However, if cycle of step size reduction exceeds 3 (see Section 
4.3.2), a new Jacobian will also be formed. 

IF JSUBITERATION is specified then program will stop after completing 
subiteration phase of solution. No ITMX card is required. 

IF SDATA CHECK is specified then program will set up network mesh 
points only. No ITMX card is required. An output file (TAPE14) with 
the network mesh points will be created for external graphics, 
processing. (Also see statement after card 22-END OF CASE). 


17 

1-10 


SPRINT 

18 

1-10 

ITPRIN 

Printing output occurs at every ITPRIN iterations. 


11-20 

ITVRCP 

= 1. for printout of variables, residuals and 
corrections resulting from full iteration 


21-30 

IPLOTP 

= 1. for printer plot of cuts of vortex system 


31-40 

IPTIME 

= 1. for printout of elapsed CPU time from 
various programs and subroutines 


41-50 

IPNPIC 

= 1. for printout of near field and far field 
information 


51-60 

IPSOLV 

= 1. for printout of out-of-core equations solver 
information 

19 

1-10 

ISINGS 

= 1. for printout of resultant values of 

singularity strength and gradient at panel 
corners, centers, and edge midpoints 


11-20 

IGEOMP 

= 1. for printout of geometry diagnostic data 


21-30 

ISINGP 

= 1. for printout of singularity spline 
diagnostic data 


31-40 

ICONTP 

= 1. for printout of control points diagnostic 
data 


41-50 

IBCONP 

= 1. printout of boundary condition diagnostic 
data (not used) 

20 

51-60 

lEDGEP 

= 1. for printout of edge matching diagnostic data 
^NETWORK 

21 

1-10 

NNETT 

Total number of networks; each call of ^VORTEX 
counts two or four networks (see ^VORTEX), NNETT 
1 20 

******************************************************************************** 


Each network is now defined in turn by a network data block which is 


headed 

by one of 

the preprocessor options 
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Data Card Variable (Numeric data are input in 6E10.0 format) 

Number Column Name Description /Comment 

^POINTS 

2QUADR I LATERAL 

^GOTHIC 

JVORTEX 

^TRAILING WAKE 

The sequence in which the networks are to be input into the program 
is irrelevant except for the following restrictions. 

1) The sequence number of the networks must be in proper 
sequential order. 

2) A network which will be updated must be input after a 
network to which it is attached. 

3) Moreover, the data cards {2QUA0R I LATERAL, ^GOTHIC, ^POINTS) 
for setting up wing and/or body networks should precede 
those (2V0RTEX, ^TRAILING WAKE, ^POINTS) for generating 
free and fed sheets and the trailing wakes. 

******************************************************************************* 
Input data for the complete case is terminated by the following card 

22 SIEND OF CASE 

If 3DATA CHECK is used, additional sets of data cases can be input 
following immediately the 3END OF CASE card. 

******************************************************************************** 

PI ^POINTS 

An input format for x,y,z coordinates of all corner points of a 
network is provided for a general cambered wing geometry or any other 
network such as body, special wake, or vortex sheets. 


H 





n 5l 9l 131 

(3) SIDE 1 ® 

N ^ 

NCOL 


CORNER NUMBERS 
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Data Card Variable (Numeric data are input in 6E10.0 format 


Number 

Column 

Name 

Description/Comment 

P2 

1-10 

KN 

Network no. 


11-20 

NT 

Type of network 


21-30 

NUP 

Update index 


= 0. Fixed network 

1. Trailing wake enamating from free or fed 
sheet 

2. Fed sheet 

3. Free sheet to which no fed sheet is attached 

4. Free sheet to which a fed sheet is attached 

P3 1-10 NROW Number of rows and columns of the specified 

11-20 NCOL network 


P4 1-60 ZM(1,I,J) x,y,z coordinates of corner points input by 

ZM(2,I,J) column (J=l, NCOL). Corner points are input 
ZM(3,I,J) sequentially (see figure above), two points per 
(1=1, NROW) card and continuous per point column. Start each 
column on a new card. 

******************************************************************************** 

Q1 1-4 $QUADRILATERAL 

This data card calls for using Quadrilateral 
preprocessor to generate mesh points for a 
specified network no. with the given four corner 
points. 
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Data Card Variable (Numeric data are input in 6E10.0 format) 
Number Column Name 


Q2 1-10 KN 

11-20 NT 

21-30 NUP 


Q3 1-60 SC(1,J) 

SC(2,J) 

SC(3, J) 

( 0 = 1 , 2 ) 

Q4 1-60 SC(1,J) 

SC(2, J) 

SC(3,J) 

(0=3,4) 

If any z coordinate of the given four corner points is not zero, a 
simple twisted wing will be set up using linear interpolation. 


Q5 

1-10 

NROW 

Number of rows 



Q6 

1-60 

YPC(I) 
(1=1, NROW) 

Percent values (100^ = !•) tor 
column, i.e., side 2 and 4. 

cuts 

along 

Q7 

1-20 

NCOL 

Number of columns 



Q8 

1-60 

XPC(J) 

(J=1,NC0L) 

Percent values (100^ = 1.) for 
row, i.e., side 1 and 3. 

cuts 

along 

Q9 



^CAMBERED WING (OPTIONAL) 




Input ^CAMBERED WING data block for simplified camber definition. If 
this option is chosen, the network mesh points generated by ?QUAD 
must have corner 1 at the apex. 


Descript ion /Comment 

Network no. 

Type of network 
Update index 

= 0. Fixed network 

1. Trailing wake enanating from free or fed 
sheet 

2. Fed Sheet 

3. Free sheet to which no fed sheet is attached 

4. Free sheet to which a fed sheet is attached 

x,y,z coordinates of corner 

pts. 1 and 2 input by the following order, 

X2*y2*^2’ 

x,y,z coordinates of corner pts. 3 and 4 
input by the following order, 
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NCEN - NCOL 


G2 1-10 KN Network no. 

11-20 NT Type of network 

21-30 NUP Update index 

= 0, Fixed network 

1. Trailing wake emanating from free or fed 
sheet 

2. Fed sheet 

3. Free sheet to which no fed sheet is attached 

4. Free sheet to which a fed sheet is attached 

G3 1-10 NCOL Number of corner points along leading edge 

G4 1-60 SC(1,J) x,y,z coordinates of corner pts. along 

SC(2,J) leading edge from nose to tail 
SC(3,J) 

(J=1,NC0L) 

G5 1-10 NROW Number of spanwise cuts 

G6 1-60 YPC(I) Percent values (100% =1.) for spanwise cuts. 

(1=1, NROW) 

G7 1-10 NCEN Number of corner pts. of wing network along 

centerline; NCEN should be less or equal to 
NCOL. If NCEN<NC0L, then wing geometry with 
swept trailing edge will result. This option is 
presently invalid for analysis due to abutment 
restrictions. It can be used to generate a set 
of data which upon proper manipulation can be 
reinput using ^POINTS. 

G8 ^CAMBERED WING (OPTIONAL) 

Input ^CAMBERED wing data block for simplified camber definition 
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Data Card Variable (Numeric data are input in 6E10.0 format) 


Number Column Name Descript ion /Comment 


T1 



^TRAILING WAKE 

The data card calls for using Trailing Wake 
preprocessor to generate mesh points for the 
trailing wake network attached to wing 




T2 

1-10 

KN 

Network no. 


11-20 

NT 

Type of network 


21-30 

NUP 

Update index 


= 0. Fixed network 

1. Trailing wake emanating from free or fed 
sheet 

2. Fed sheet 

3. Free sheet to which no fed sheet is attached 

4. Free sheet to which a fed sheet it attached 


T3 

1-10 

NAT 

Sequence number of network to which side 1 of 
trailing network is attached 


11-20 

NSD 

Side of network to which side 1 of trailing 
network is attached 

T4 

1-10 

XWAKE 

X coordinate of corner pt. of downstream of the 
trailing wake; should be about 50 times of the 
X-coordinates of the trailing edge along 
centerline. It is essential that this value 
should be the same as the one given in card no. 


V6 under 3V0RTEX. 

******************************************************************************* 
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SIDE 3 



Data Card Variable (Numeric data are input in 6E10.0 format) 


Number Column Name Descript ion /Comment 


VI ? VORT EX 

This data card calls for using Vortex 
preprocessing to generate mesh points for free 
sheet, fed sheet, and as an option the attached 
trailing wakes. Either two or four networks will 
be generated. 



V2 

1-10 

KN 

Network no. for 

free sheet, the fed sheet will 




have network no. 

KN+1. 


11-20 

KW 

If KW = 1., then 

the trailing wake networks 


attached to free sheet and fed sheet will also be 
formed. Their network numbers will be KN+2 and 
KN+3 respectively. 

If KW=0. then no trailing wake networks are formed 


V3 

1-10 

NROWF 

Number of rows on free sheet 




11-20 

NROWV 

Number of rows on fed sheet 



V4 

1-10 

NAT 

Sequence number or network to 
the free sheet is attached 

wh i ch s i de 1 

of 


11-20 

NSD 

Side of network which side 1 

of free sheet 

is 


attached 
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Data Card Variable (Numeric data are input in 6E10.0 format) 


Number Column Name 


Descript ion /Comment 


V5 


V6 


21-30 NAIF 

I- 10 ARC 

II- 20 JNAT 

1-10 XWAKE 


Sequence number of another free sheet network to 
which side 4 of the free sheet is attached 
Set NATF=0 if side collapses to a point 

Perturbation parameters for size of free and fed 
sheets. The parameter a is reset to a + ARC. If 
the user wants the original initial guess, ARC 
should be set to 0. 

Sequence number of the vortex network to which 
the apex of the complete vortex system resides. 
JNAT will differ from KN when more than one set 
of networks are used to define the free and fed 
sheets. See figures 25 and 26 as examples. 

This card is required only if KW = 1. 

X-coordi nates of corner pt. at downstream of the 
trailng wake; should be about 50 times of the 
X-coordi nate of the trailing edge along 
centerline, (also see T4) 


******************************************************************************* 


^CAMBERED WING 

A deck must first be prepared to generate the desired networks for the flat 
plate representation of the configuration to be studied. The wing plan view 
itself will be generated either through use of the ^QUADRILATERAL preprocessor 
or the 2G0THIC preprocessor. 

The three-dimensional character of the wing can be defined by use of the 
$CAMBERED WING preprocessor. This preprocessor generates the z coordinate for 
the (x,y) coordinates of the flat wing representation of the desired 3-D wing 
through interpolation. In general, the cambered surface is defined through a 
set of input data specifying the wing mean lines in the chordwise direction at 
a limited number of spanwise stations (no more than 50). It is also possible 
to input a fixed mean line shape valid for all span stations scaled to the 
local chord. The 3CAMBERED WING preprocessor also can generate the camber 
surface for wings with circular arc spanwise camber. This preprocessor was 
originally developed in reference 6. 

Thus, the current technique for generation for three-dimensional wing networks 
consists of two steps: 

(1) Generate wing plan view, with desired paneling density using ^GOTHIC 
or 2QUADR I LATERAL. 

(2) Generate wing z coordinates using 55CAMBERED WING. 
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It is essential that the card set $CAMBERED WING should follow immediately the 
card set $QUADRIALATERAL or $G0THIC. 

A description of the input card preparation, as part of the $G0THIC or 
$QUADR I LATERAL input cards is as follows (data are input in 6E10.0 format): 


Number Column Name 


Descri ption/Comment 


Cl 


$CAMBERED WING 


C2 1-10 CNTRL CNTRL controls which type of wing is 

generated: CNTRL = 1. is for a single mean 

line for all span stations, CNTRL =2. is for 
varying camber and twist with span, CNTRL = 3. 
generates the Wentz (ref. 5) conical cambered 
delta, and CNTRL = 4. generates the Barsby 
conical cambered deltas (see equation (15)). 


The input cards hereafter differ and will be described for each of the 4 
possible values of CNTRL. 


(1) If CNTRL = 1, a 3-0 wing with a single camber shape will be 
generated. The necessary input cards are as follows: 


C3 

C4 


C5 


1-10 

NPCT 

Number of x/c's (of which z/c's will be 
defined) on card C4. 

1-60 

PCTX(I) 

PCTX(I) is a table of percent local 


(1=1, NPCT) 

chord at which the z percent local chord is to 
be specified on the following cards. (100% = 
1.0) 

1-60 

PCTZ(I) 
(1=1, NPCT) 

PCTZ(I) is a table of the z values in z/c. 


This completes the necessary input for a general wing with a single mean 
line shape. The desired z values for the paneling generated in SGOTHIC or 
$QUADR I LATERAL are then found through linear interpolation. 


(2) If CNTRL = 2., a 3-D wing with camber and twist varying with span 
station will be generated. (Note that for CNTRL = 1. or 2., it is 
not necessary to specify NPCT = NROW or NTST = NCOL for the network 
in question.) The only restrictions are NPCT, NYST > 50. The 
necessary input cards are: - 


Number 

Column 

Name 

Description/Comment 


C3 

1-10 

NPCTl 

Number of x/c's at which 
card C4 

z/c's will be defined on 


11-20 

NYST 

Number of y stations at 
defined, the x/c array 
station 

which z/c's will be 
will apply at each y 



Number 

Column 

Name 

Descr i pt i on/Conment 

C4 

1-60 

PCTX(I) 

(I-1,NPCT) 

PCTX(I) is a table of percent local 
chord at which the z percent local chord is 
to be specified on the following cards. 
(100% = 1.0) 

C5 

1-10 

YSTA 

Y-location at which array of z/c's will be 
defined. Input NYST sets of C5 and C6 data 
cards 

C6 

1-60 

PCTZ(I) 

(I=1,NPCT) 

PCTZ(I) is a table of the z values 
in z/c. 


This completes the necessary input for a general wing with varying twist 
and camber shape. The desired z values for the paneling generated in 
$G0THIC or $QUADR I LATERAL are then found through linear interpolation. 


(3) If CNTRL = 3., a conically cambered delta wing will be generated 
where the first (0.805) b/2 is flat and at the maximum z, and the 
remainder of the wing semi span is a portion of a circular arc. The 
maximum z is 0.105 of the wing local semi span. See reference 5 for a 
description of these wings. When CNTRL = 3., jto further data cards 
are required. 


(4) If CNTRL = 4., a conically cambered delta wing will be generated, 
where the wing is a portion of a circular arc in the spanwise 
direction determined by the equation: 


(15) 


’1 ocal 


(b/2) 


local 


(1 + P^) 


2p y 


JW2j 


local 


local 


-yo + 


p') - (2P)^ 


where p = 0.0 corresponds to a flat wing an p = 1.0 corresponds to a 
wing which is one-half of a cone. One further card is then required 
to specify the value of p. 


Number Column Name 


Descri pti on/Comment 


C3 1-10 pi 


0.0< pl< 1.0 
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5.6 Example Input Case 


As an aid to the user in understanding the proper application of the input 
format specifications, two example input cases are provided. The first 
case consists of an aspect ratio 1.15 flat delta wing with a design wake. 
The network arrangement is shown in Figure 25, section 5.4.2. The inputs 
for this case are given in Figure 31. 

The second case is a 70 panel twisted arrow wing. The network arrangement 
is shown in Figure 26, section 5.4.3. The inputs for this case are given 
in Figure 32. 

5.7. Practical Instructions 


The preceding sections have given sufficient instructions to properly set 

up a solution, here practical hints are given to aid in their use. 

1. Always submit a data check (JDATA CHECK, Card 15) before submitting 
the solution to iteration. Check abutment data (described in section 
6.1.3) to ensure proper network arrangement. 

2. Use ^SUBITERATION (Card 15) for further checking of unusual cases. 

If solution does not converge in subiteration it will not converge in 
the full iteration process. Use of JSUBITERATION is not necessary 
for typical cases. 

3. Most well posed cases using conical type paneling on the lifting 
surface will converge using ^ITERATION (Card 15). Try 5 iterations 
(ITMX, Card 16) and save results on TAPE14, If case appears to be 
converging but residual (SSR) are greater than 10“^, repeat 
iterations starting with saved results (set ITMX negative). 

4. Cases which appear to have difficulty in converging using ^ITERATION 
may respond to the least squares solver, 2LEAST SQUARES. This will 
include most cases with streamwise lifting systems. Save data on 
TAPE14 in case more iterations are necessary. Remember that JLEAST 
SQUARES is more expensive than ^ITERATION. Also remember the 
paneling limitations given in Table 1. 

5. Reasonable results have been obtained using 60-70 panels on the wing, 
7-9 rows (NROWF, Card V3) on the free sheet, and 3 rows (NROWV, Card 
V3) on the fed sheet. More wing panels may be necessary to obtain 
desired simulation of wing camber and resolutuion of the pressure 
distribution. 
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FIGURE 31 INPUT SPECIFICATION - AR = 1.15 DELTA WING 



A COMPUTER PROGRAM 
FOR 

A THREE DIMENSIONAL SOLUTION OF FLOWS OVER WINGS 
WITH LEADING EDGE VORTEX SEPARATION 


- LIST OF INPUT data CARDS - 


NO. CARD images 

1 SCASE 

2 70 PANEL TWISTED ARROW WING OF MANROE 

3 TINOCO AND LU 

A JANGLE OF ATTACK 

5 15.9 

fc SYAW 

7 0. 

e jeymmetry 

9 1. 


1 0 

JMACH NO. 






1 1 

0.40 






] 2 

JRFFERENCES 






13 

26.777 0. 

0. 





14 

4C9.713 1. 

29 

.652 1 

.0 



15 

tDATA CHECK 






16 

JPRINT 






17 

1. 1. 

1. 





1 A 

0. 0. 

0. 

0 

• 6 * 

0. 


19 

SNFTWnpK 






2 : 

1C. 






21 

JPOINTS 






2 2 

1. 2. 

0. 





23 

8. 






24 

3.51020 

0.00000 

0.00000 

3.51020 

0.00000 

0.00000 

25 

3.51020 

o.oocco 

0.00000 

3.51U20 

0.00000 

0.0 JOOO 


3.51020 

0 . 0 G C C 0 

0.0 COCO 

3.51 ..20 

0 . G C u G 0 

0. C 0000 

27 

3.51020 

O.GOCOO 

0.00000 

3.51020 

0.0 0 000 

0.00000 

26 

9.76610 

0.00000 

o.ocooo 

9.76610 

.29813 

-.00992 

29 

9.7661C 

.76662 

-.026 18 

9.76610 

1.21382 

-. 04248 

30 

9.76610 

1.51195 

-.05377 

9.76610 

1.72490 

-. 06204 

31 

9.76610 

1.95914 

-.07133 

9.76610 

2.12950 

-. 09286 

32 

la. 00000 

0.00000 

O.GOOOO 

18.00000 

.69054 

-.01605 

33 

16.00000 

1.77566 

-.04493 

18.00000 

2.81147 

-. 14472 

34 

18.00000 

3.50200 

-.23210 

18.00000 

3.99524 

-.29902 

35 

18.00000 

4.53781 

-.3E9C5 

18.00000 

4.93240 

-.45796 

36 

26.E00CG 

o.ooooc 

O.OCOOO 

26.6GCC0 

1.10039 

-.C1407 

37 

26.60000 

2.82956 

-.06644 

26.60000 

4.48014 

-.23916 

38 

26.60000 

5.58053 

-.37778 

26.60000 

6.36652 

-.49201 

39 

26.60000 

7.23111 

-.63322 

26.60000 

7.85990 

-.74483 

40 

35.00000 

0.00000 

O.OUOOO 

35.UO00D 

1.50070 

-.00385 

41 

35.00000 

3.85895 

-.06739 

35.00000 

6.11000 

-.22973 

42 

35.00000 

7.61010 

-.39398 

35.00000 

8.68263 

-.53159 


FIGURE 32 INPUT SPECIFICATION - ARROW WING 
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47 

35.00000 

9.86176 -.70224 

35.00000 

10.71930 


-.82589 

44 

40.00000 

0.00000 0.00000 

40.00000 

1.73898 


. 00608 

45 

40.000CD 

4,47167 -.02067 

40.00000 

7.08014 


-.16643 

46 

40.00000 

8,81912 -.33143 

40 .00000 

10.06125 


-.47759 

47 

40.00000 

11,42760 -.64781 

40.00000 

12.42130 


-.77641 

48 

45.21610 

0.00000 0,00000 

45.21610 

1.98757 


.01949 

49 

45.21610 

5,11066 .06109 

45.21610 

8.09223 


-. 05161 

50 

45.21610 

10,07960 -,21367 

45.21610 

11.49949 


-.35829 

51 

45.21610 

13.06115 -.53012 

45.21610 

14.19690 


-.65902 

52 

»POINTS 






53 

2 . 2 • 

0. 





54 

8. 5. 






55 

45.21610 

0.00 COO o.ncooo 

45.21610 

1.98757 


. 01949 

56 

45.21 610 

5.J10P6 .06109 

45.21610 

0.09223 


-. 05161 

= 7 

45.21610 

10.Q7980 -.21367 

45.21610 

1 1.49949 


-.35829 


45.21 610 

13.06115 -.53012 

45.21610 

14.19690 


-.65902 

r o 

50.00000 

6.01431 .UfcOl 

50.00000 

6.52785 


,16834 

60 

5O.?0Q&C 

8.90626 .10163 

50.00000 

11.17657 


-• 05564 

6 1 

50.CP0 12.69011 -.19566 

50.00 

13.77121 


-.30487 

6? 

50.0(i000 

14,96042 -.42570 

50.00000 

1 5.02530 


-.51511 

63 

55.00000 

10,25512 .27657 

55.00000 

11.27323 


, 21775 

64 

55.00000 

12.87311 .09550 

55. ncooo 

14,4 0026 


-.04357 

65 

55, CO C JO 

15.41837 -.14196 

55.00000 

16,14559 


-.21413 

6 6 

55.C0C.0O 

16.9455* -.29509 

55.00000 

1 7.52730 


-.35508 

6 7 

60. OOOOO 

15.4959-* .18161 

60.00000 

1 6.01062 


.13249 

6b 

6D.OOOOO 

16.83999 .05359 

60.00000 

17,62401 


-. 02349 

69 

6Q.OC000 

18.14670 -.07615 

60.00600 

10,52004 


-.11541 

7 n 

6'.ocoor 

18.93072 -.15998 

60,00000 

19,22940 


-.19332 

71 

05.27500 

21.02500 -.00004 

65.27500 

21.025C0 


-.0'J304 

72 

65.27500 

21.02500 -.00004 

65.27500 

21.02500 


-.00004 

73 

65,275' 0 

21.02500 -.00004 

65.27500 

21.02500 


-.00004 

74 

65.275uO 

21.02600 -.00004 

65.27500 

21.02500 


-.00004 

7 5 

tVORTEX 






76 

3. 0. 
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8. 3. 






7« 

1 . 3 . 

0. 
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0. 3. 






PO 

IVORTEX 
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B3 

2. 3. 

3. 





B4 

0. 3. 
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2500. 
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87 

9. 6. 

0. 





88 

4 . 5 . 






09 

45.2161 0. 

0. 

52.5 

0. 

0 

• 

90 

6 0, 0 • 

0. 

65.275 

0. 

0 

« 

91 

50. 5. 

01431 .15601 

55. 

5.1014 


156 

92 

6U , 5 . 

014 .156 

65.275 

5.C14 


156 

93 

55. 10 

.25512 .27657 

57.5 

10.26 


277 

94 

62,5 10 

.26 .277 

65.275 

10.26 


277 

95 

60. 15 

.49594 .18161 

62. 

15.50 


182 

96 

64 . IS 

.50 .182 

65.275 

15.50 


182 

97 

65.275 21 

.025 -.00004 

65.275 

21. 025 

- 

.00004 

9fc 

65,275 21 

.025 -.00004 

65.275 

21.025 

- 

.00004 

99 
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100 

10. 8. 

0. 
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9, 3. 






ii.2 

2500. 






103 

lEND OF CASE 






FIGURE 32 CONCLUSION 


6.0 OUTPUT GUIDE 


In this section the organization of the computer output is described. The 
nomenclature employed in the output is summarized in Table 4. A typical 
output (with the appropriate print options) will consist of a copy of the 
program inputs, printer plots of the initial free/fed sheet, network 
coordinates, and network abutment data. If the JDATA CHECK option had been 
used the program would terminate at this point. If the JSUBITERATION option 
is used then printout will also include results from the initial solution. If 
a complete solution is sought the singularity parameters, corrections, and 
residuals, and updated geometry can be printed every iteration. Printout of 
the solution results (pressure coefficients, velocities, etc) can be deferred 
until the final iteration, printed every so many iterations, or every 
iteration. 

Additional print options provide for the printing of detailed diagnostics 
which can be of use to a user intimately familiar with the workings of the 
code. Although of not much use to the typical user a brief description will 
be given of these options. 

The program also creates a file (TAPE14) which can be saved and used to 

restart the iterations if convergence is not achieved in the first solution 

attempt. A description will be given of TAPE14 in section 6.4. 

Examples of the output will be presented to aid in its description. The 

examples used are from the aspect ratio 1.15 flat delta wing case used to 
illustrate the inputs in section 5.6. Since the program always prints out the 
inputs, the first part of the output will appear as shown in Figure 31. The 
input formats are described in section 5.5. 

6.1 Data Check Output 


In this section the print out typically associated with a data check will 
be described. This print out is included in every solution. 

6.1.1 Free and Fed Sheet Printer Plots (OPTIONAL, IPL0TP=1) 


Printer plots of the free and fed sheets geometry are shown in Figure 33. 
These plots are created when IPL0TP=1 (Card 18). The network arrangement. 
Figure 25, contained two sets of free and fed sheet networks which results 
in the two plots. A different symbol is used for each transverse cut. 

The poor resolution of the printer plots results in a somewhat ragged look 
of the cuts but will at least give the user an idea of the initial 
geometry. Note the change in the horizontal scale for the plots of 
networks 6 and 7. 

6.1.2 Mesh Point Data 

The program always prints the mesh point data. An example of the mesh 
point data printout is shown in Figure 34. Only a partial list is shown 
for the sake of brevity. The mesh point data is organized by network. 
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INPUT nCTudrk ntsh points data 



1 

IHUMBCR pout 

s 7 NUMBER COlUMNI b 61 

••ODOCt 

O.COObO 

e.oceoo 

0.00000 

0.00000 

O.OOOOO 

1. 00(00 

c.sceco 

D.ocooe 

0 . ocooo 

0.00000 

o.ooooo 

i.oosoe 

0.00000 

o.oooco 

o.oeeoo 

0.00000 

0.00000 

o.eoMt 

0.00000 

0.00000 

.20000 

o.ooooo 

0.00000 

• ZOOM 

.01726 

0.00000 

.20000 

.82875 

e. 00000 

.20001 

.03544 

0.00000 

.20000 

•04313 

0.00000 

•2000V 

. 9SD31 

0.00000 

.20000 

.05750 

o.oocoo 

.40000 

o.tuooc 

0.0C09C 

.4C000 

.93450 

o.oocoo 

.40000 

.09750 

o.oooco 

.40000 

.07188 

0.00000 

.4CC00 

.08625 

9.GC030 

.40060 

.16063 

o.coooo 

.400CO 

.11500' 

0.03000 

.60090 

o.cQfoe 

0. 00000 

.SDOtPO 

.08175 

0.00000 

.60060 

.08625 

0.00000 

.60000 

.10781 

0.00000 

.60000 

.12938 

0.00000 

.60000 

.15094 

o.oocoo 

.60000 

.17290 

0.00000 

.80000 

0. 00000 

0.00000 

.80000 

.06900 

0.00100 

.00000 

.11500 

0.00000 

.00000 

.14375 

0.00000 

.80000 

.17250 

0.09000 

.6oa"o 

.20129 

o.oocoo 

.80000 

.23000 

0.00000 

1.00000 

O.OOOCO 

0.00000 

1.00000 

.08629 

0.00000 

1. 80000 

.14375 

0.00600 

1.00000 

.17969 

o.oocoo 

1.00000 

.21563 

o.rocoe 

l.OOODt 

.25156 

0.00000 

1 . ocooo 

.20750 

0.00060 


NETWORK NO. 

2 

(NUMBER ROWt 

= 3 1 

NUMBER COLUMNS * 7| 

l.OPOOl 

0. OCOOP 

0 . ocnco 

1. 65000 

O.OOOOO 

0.00000 

1.10006 

0. 00000 

C. 00000 

1.00000 

.08625 

0.00060 

1 .05000 

.09956 

D.OPPOC 

1.10600 

. 0948 R 

0.00060 

1 .60000 

.14375 

0 .POO nc 

1 .05CP0 

.15094 

D.ggr -g 

1.10000 

.15013 

o.oooco 

i.oDooe 

.17969 

0.00600 

i.oSooo 

. 18817 

O.OOOCO 

1. 10000 

.19766 

1.00000 

l.OOOCI 

.21563 

0.00000 

I.OSCOO 

.22641 

o.oooer 

1.10000 

.23719 

0.00000 

i.Dooeo 

.25156 

o.ooooo 

1 .05000 

.26414 

0.09000 

1.10000 

.27672 

o.ooooo 

i.ooboe 

.28750 

O.C9CCO 

1.05090 

.36108 

0.00000 

1.10000 

.31625 

0.06000 





HETUORR NO. 

3 

INUBSCP ROW! 

= 2 

NUMBER COLUMNS s 7| 

1.10000 

0.00000 

0.00000 

50.00000 

O.OOOOO 

0.00000 

1.10600 

.09480 

0.69600 

50.00000 

.094(0 

o.oooco 

1.10006 

.15813 

0.00000 

sc.ooonc 

.15813 

0.00000 

1.10600 

.19766 

0.00000 

50.00000 

.19766 

0.00000 

i.toooo 

.23719 

0.00000 

50.00000 

.23719 

o.ooooo 

i.iDooe 

.27872 

O.OOOCO 

50.06000 

.27672 

0.00000 

i.iDoeo 

.31625 

o.ooooo 

SO.OOODO 

.31625 

> 0.00000 


NETWORK NO. 


(NUM8CR ROUS 

c 9 NUMBER COLUMNS c 61 

0.06060 

0 . 0 (1 0 0 c 

0.66000 

0.06000 

0 . ocooo 

0.00008 

o.occto 

0.9CP0O 

O.OOOCO 

0.00000 

0.00000 

0.06000 

0.00060 

0.00000 

O.DOOCQ 

0.00000 

0.00000 

O.OOOOO 

O.OOOOO 

0.00000 

0.60000 

0.00000 

0.00000 

o.ooooo 

o.ooooo 

0.9(000 

0. 00000 

.20000 

.05758 

0 . oooco 

.20000 

.06211 

.00563 

.20000 

.08398 

.01277 

.20000 

. 08392 

•62015 

.20000 

.06241 

.02731 

.20000 

.05969 

.034K 

.20090 

.05545 

.64015 

.20000 

.04943 

•04439 

*20000 

.04233 

.04617 

.4CC00 

.11500 

0. 00(00 

.40000 

.12472 

.01126 

.40000 

. 12796 

•02554 

.40000 

.12783 

.04030 

.40000 

.12483 

.05461 

.40000 

.11938 

.06820 

•40000 

.11091 

. 68029 

.40000 

.09685 

.08677 

.40000 

.00467 

.09233 

.60000 

.17250 

G.OUbCP 

.60000 

.1B632 

.dlC89 

.60000 

.19195 

.03831 

.60000 

.19175 

.06845 

.60000 

.18724 

.68192 

.CP6CI 

.17906 

.10230 

.60000 

.18635 

.12044 

.60000 

.14828 

.13316 

.60000 

.12710 

.13850 

.80000 

.23600 

0.09000 

.86000 

.24843 

.02251 

.80000 

.25593 

•05108 

.80000 

.25567 

.08061 

•00000 

.24965 

.-10923 

.00000 

.23075 

.13640 


FIGURE 34 MESH POINT DATA 



Each network is identified followed by the number of rows and columns (M 
and N), Paired mesh points coordinates are presented by row in 6F10.5 
format in X Y Z order. 

6.1.3 Abutment Data 


The program always prints the abutment list, abutment intersection list 
and the update index arrays. These data are vital to determining whether 
the networks have been properly defined to ensure the appropriate matching 
along their edges. Unless proper matching occurs along all network edges 
the flow model and the resulting solution will be in error. Every network 
arrangement should be run through data check and have its abutment list 
thoroughly checked before being committed to solution. 

The abutment list is shown in its entirety in Figure 35. The 
corresponding network arrangement is illustrated in Figure 36. A similar 
sketch should always be prepared to aid the user in the abutment checks. 
The abutment list consists of the abutment number, side and network 
involved, and a characterization of the control points. The user should 
refer to section 5.2.4 for the control point placement for the various 
types of networks. To facilitate understanding of the abutment list a 
walk through of the list follows. 

Abutment 1 concerns side 1 of network 1 which is on the plane of 
symmetry. These control points satisfy their appropriate boundary 
conditions. Side 4 of network 1 collapses to a point and has no 
abutment. Abutments 2 and 3 concern the remaining two sides of network 
1. These sides both abut to other networks, side 2 to the near wake 
network 2 and side 3 to the free sheet network 4. Proper matching is 
indicated by the characterization "control points peform doublet 
matching." In both cases the control points on network 1 satisy the 
boundary conditions while the control points on networks 2 and 4 do the 
matching. Abutment 4 between networks 2 and 6 is similar. 

Abutments 5, 10, 11, 14, and 15 are between networks in which one edge has 
no control points. In each case the proper abutment is indicated by the 
edge with the control points performing doublet matching. Abutments 7, 

12, 16 and 18 have no control points on their coincident edges. Here 
proper abutments are indicated by the network pairing. Abutments 8, 13, 
17, 19, 20 and 21 are not actually abutments but free edges. If an error 
is made in the network definition such that the proper abutment is not 
made between two adjacent networks, the program will not pair the networks 
and regard them as free edges. 

The Abutment Intersection List is shown in Figure 37. This list is 
similar to the abutment list except that it characterizes the control 
point behavior at network corners. If the abutment list checks out for 
the configuration the user need not concern himself with the intersection 
list. 

The Update Index Arrays are shown in Figure 38. These arrays list the 
various control parameters which regulate the geometry update during 
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AtUTPEHT IIST 


SIOI 


CHMtCUItZATION 


Control points iif anti use oriciml rouudart conditions 


CONTROL POINTS UP ANTI USE 0RI6INAL OOUNDART CONOITIONS 
CONTROL POINTS PLRFORN OOURLET HATCHIN6 


CONTROL POINTS (IP ANTI USE OftlOlNAL NOUADAOT CONDITIONS 
CONTROL POINTS PERFORH OOUOLET NATCHINS 


CONTROL POINTS (IF ANT) USE ORIOINAL iOUNOART CONDITIONS 
CONTROL POINTS PERFORH OOUOLET NATCHINS 


NO CONTROL POINTS 

CONTROL POINTS PCNFORN DOURLCT NATCHINS 


CONTROL POINTS (IF ANTI USE ORIStNAL OOUNOART CONDITIONS 


NO CONTROL POINTS 
NO CONTROL POINTS 


NO CONTROL POINTS 


CONTROL POINTS IIF ANT) USC ORISINAL SOUNDART CONOITIONS 


NO CONTROL POINTS 

CONTROL POINTS PERFORH OOUOLET NATCHINS 


NO CONTROL POINTS 

CONTROL POINTS PERFORH DOUBLET NATCHINS 


HO CONTROL POINTS 
NO CONTROL POINTS 


NO CONTROL POINTS 


NO CONTROL POINTS 

CONTROL POINTS PERFORH DOUBLET. NATCHINS 


NO CONTROL POINTS 

CONTROL POINTS PERFORM DOUBLET NATCHINS 


HO 

CONTROL 

POINTS 

NO 

CONTROL 

POINTS 

NO 

CONTROL 

POINTS 

NO 

CONTROL 

POINTS 

NO 

CONTROL 

POINTS 

NO 

CONTROL 

POINTS 


FIGURE 35 ABUTMENT LIST PRINTOUT 
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COMTROL POINT PERFORHS DOUBLET H4TCHIN0 ON SIDE 3 


CONTROL POINTS (IF ANTI USE ORIOINAL lOUNDART CONDITIONS 
CONTROL POINT PERFORHS DOUBLET HATCHINt ON SIDE 1 


CONTROL POINT PERFORHS DOUBLET NATCHINO ON SIDE 3 
CONTROL POINT PERFORHS DOUBLET HATCHINO ON SIDE I 
NO CONTROL POINT 

CONTROL POINT PERFORHS DOUBLET HATCHINO ON SIDE I 


NO Control point 

CONTROL POINT PERFORHS DOUBLET PATCHINC ON SIDE 1 
NO CONTROL POINT 

CONTROL POINT PERFORHS DOUBLET HATCHINO ON BIDE 4 


NO CONTROL POINT 

CONTROL PCINT PERFORHS DOUBLET HATCHINO OH SIDE 1 


NO CONTROL POINT 

CONTROL POINT PERFORHS DOUBLET HATCHING ON SIDE 1 
NO CONTROL POINT 

CONTROL POINT PERFORHS DOUBLET HATCHING OH SIDE 4 


HO CONTROL POINT 

CONTROL POINT PERFORHS DOUBLET HATCHING ON SIDE 1 
CONTROL PtINT PERFORHS DOUBLET HATCHING ON SIDE 2 
CONTROL POINT PERFORHS DOUBLET HATCHING ON SIDE I 


FIGURE 37 ABUTMENT INTERSECTION LIST PRINTOUT 


77 


CO 


UPDATE INDEX ARRAYS 
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91782222 AlAA 7 8 1 3 1 2 3 3 9 9 


FIGURE 38 UPDATE INDEX ARRAYS PRINTOUT 



iteration. These data define how the various networks are connected 
together. Table 4 lists the definitions of the various headings. The 
general user need not be concerned with this print out. 

6.2 Solution Output 


In this section the printout typically associated with an iteration 
solution will be described. For every main iteration the following is 
printed . 

Iteration Summary 
Iteration No. 

Sum of Squares of Residuals (SSR) = 

No. of Function called = 

Fraction of Newton Step taken = 

Step Size (Length of Correction Vector) = 

Force and Moment Data for NT = 1 or NT = 2 Networks 
Network Mesh Point Data 

For every ITPRINT (Card 18) iterations the detailed physical quantities 
are printed. 

When IPLOTP = 1 (Card 18) 

Printer Plot of Vortex Sheet 
When ITVRCP = 1 (Card 18) 

Values of Variables 

Residuals 

Corrections 

6.2.1 Iterative Results Summary 

These data are printed every main iteration and serve to summarize the 
progress of the solution during the iteration process. 

Iteration No. - Counter on the number of iterations taken. Starts with 
"0" for the initial guess solution. 

Sum of Squares of Residuals (SSR) - This is the sum of the squares of the 
residuals. When using ITFLOW, SSR = F^ + (refer to equation 9, 
section 4.3.1). In using the Quasi-Newton, ITFLOW, SSR is a reliable 
indicator of the goodness of the solution. A solution may be considered 
acceptable when SSR < 10“^. 

When using the least squares solution procedure LSFLOW, SSR = 
(subiteration drives F to zero, also, refer to equation 12, section 
4.3.2). A solution is acceptable when SSR < 10“^. For higher values of 
SSR the residual is not a reliable indicator. Many converged solutions 
will have residuals larger than 10“^. This may result because the panel 
twist residual K is included in the sum. A more reliable indicator in 
this case may be the residuals of G. 
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No. of Function called - Cumulative number of function residual 
canculations for the iterations performed. 

Fraction of Newton Step Taken - Refer to Appendix G, Volume I - Theory 
Document. 


x(i+l) = 0<6<1 

where 6 is the fraction of Newton step taken and is the correction 

vector. 

Step Size (Length of Correction Vector) - Refer to Appendix G, Volume I - 
Theory Document 

|l6Ax'’'|| = 6|jAX<’^|| 
where || |j is the euclidean length. 

Force and Moment Data for NT = 1 or NT = 2 Networks - This is a force and 
moment summary on all source /an a lysis (NT = 1) and doublet/analysis (NT = 
2) type networks. These type networks are used to define conf iguration 
surfaces and this summary will give the total forces acting on the 
configuration. An example of this printout is shown in Figure 39. The 
force and moments are with respect to the configuration axis system. Also 
printed is the network surface area. Three sets of data are printed for 
each quantity. The top values represent forces and moments calculated 
integrating pressures (based on the isentropic formula, section 6.2.2, 
equation 19) on the upper side of the network, middle values represent the 
lower surface totals, and the bottom values represent the sum of upper and 
lower. The upper and lower sense is determined by the NxM vector (points 
out from the upper surface). 

Network Mesh Point Data - Updated mesh point geometry. See 6.1.2 for 
details. 

6.2.2 Detailed Physical Quantities 


For every ITPRINT (Card 18) iterations, the detailed physical quantities 
are printed. These are also printed for the initial solution and the 
final iteration. An example of these results are shown in Figure 40. The 
quantity headings are defined in Table 4. The printout is organized by 
network. For every panel center-control point the following quantities 
are listed: 

Source strength 
Doublet strength 
Doublet strength gradient 
Perturbation velocity potential* 

Total velocity potential* 

Perturbation mass flux vector* 

Total mass flux vector* 
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40 PHYSICAL QUANTITIES PRINTOUT 



Normal and tangential components of total 
Mass flux vector* 

Normal component of perturbation mass flux vector* 
Pressure coefficient* 


Quantities starred are listed for both the upper and lower surfaces of the 
network. Four different pressure formulas are used for the pressure 
coefficient: 


C 


P = 


C 


P = 


C 


P = 


c 


p = 


- 2u (Linearized) 

- 2u - v2 - (Slender Body) 

- 2u - (1 - M^) - v^ - w^ (Second Order) 


2 

' ' 'oo 





( Isentropic) 


( 16 ) 

( 17 ) 

(18) 

(19) 


Here (u,v,w) is the perturbation velocity vector, referred to the 
compressibility axis which is aligned with the freestream vector. Also 
listed are the differences in the four pressure coefficients across the 
network (upper surface value minus lower), the control point label, the 
panel label, the control point coordinates, and the components of the 
freestream velocity. The results for each network are followed by 
summaries of the forces and moments similar to those shown in Figure 39 
and described in the previous section. In addition to the total forces 
for the network, forces and moments are also given for each column of the 
network . 


6.2.3 Free and Fed Sheet Printer Plots (OPTIONAL, IPLQTP= 1) 


When IPL0TP=1 (Card 18) printer plots of the updated free and fed sheet 
are produced. These plots have been illustrated in Figure 33 and are 
discussed in section 6.1.1. 

6.2.4 Variables, Residuals, and Corrections (OPTIONAL, ITVRCP=1) 

When ITVRCP=1 (Card 18) values of the program variables, the residuals at 
each control point, and corrections for the iteration are listed. These 
data are organized by type and by network. Figure 41 illustrates a 
partial printout of these quantities. The singularity strengths, the 
corrections to these strengths, and the associated residuals are ordered 
by control point order. See section 5.2.4 on network types and uses for 
proper ordering. The orientation angles, corrections and associated 
residuals are ordered one per panel and are given in degrees. Geometry 
parameters lamda and nu, corrections, and residuals are ordered by 
free/fed sheet column. Note that the residuals are identified as to type 
and that the sum of the squares of the residuals is given for each network. 
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..T9a9C-B2 

.24R2C-R2 

.3BR4C-R2 

-.5RB7C-ai 

-.4R21C-B2 

.1391C-R2 

-.4SA9C-B2 

-•fcR2BC'-B2. 

.9IC4C-B3 
.239TE-a2 - 

-.9349C-at - 

-.4Riac-a2 - 

•t379C-B2 

.9a79C-B2 

.T192C-R2 

.2«17C*a2 

.442RE-a2 

-.Asasc-at 

-.«724C-at 

■1919C-R2 

-.9a31C-B2 

•.4RR3C-RC 

••244RC-R2 

.Il«9C-a2 

-•3a22C-B2 

-..C9S4C>BI 

.33aac-af 

>.3aisc-a.a 

INFCRRCABLC R.C. FOR HI 

ITHORK NO. 

4 


tasH 

« .3792C' 

-aai 



-.ia99C-at 

.4794C-82 

-.223CC-R1 

-.5399C-at 

-.1994E*at 

-.3Raac-Ri 

-.1449C-R1 

•9aR9E-R2 

•.a99HE-R2 

.U4BC-R1 

.4RT3C-B2 

.TR4tC-a2 

.S124C-R2 

.9947C-at 

•.1443C-R1 

-.25aic-ai 

FORCE R.C. 


NO. 9 



IRRR 

« .T949C 

>031 



.1B79C-B1 

r.49aiC-R2 

• n93C*Bl 

-.4499C>-a2 

..iRRac-ai 

-.39R2C-R2 

• liRac-Ri 

-•3421C-R2 

• 1244C-R1 

-•oaotc-oi 


FORCE R.C. FOR HCTHORK NO. 7 ItRR > .IROaEvRII 

>.13tlC-«2 ••3R74C-H •’•a771C-;tt •.333(C'«« 

FIGURE 41 CONCLUDED 
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6.3 Diagnostic Printouts 


In this section the various optional diagnostic printouts will be briefly 
described. These printouts were used in the early program development and 
are as such not much use to the typical user. 

6.3.1 Geometry Data (OPTIONAL, I6EQMP=1) 


An alternate form of the geometry data is printed when IGE0MP=1, (Card 
19). An example of this printout is shown in Figure 42. The program 
lists the coordinates of all grid points along with its number, row 
number, column number, and network number. Panel data is also included 
listing panel number. 

6.3.2 Singularity Distribution Definitions (OPTIONAL, ISINGP=1) 

Panel distribution quantities shown in Figure 43 are printed when ISINGP=1 
(Card 19). This printout lists the singularity parameter numbers and 
coefficients for the nine canonical points on each panel. Refer to 
Appendix B, Volume I - Theory Document, for further discussion of these 
parameters. 

6.3.3 Control Point Data (OPTIONAL, IC0NTP=1) 


Control point data, shown in Figure 44, are printed when IC0NTP=1 (Card 
19). These data list the control point index, network number of control 
point, panel number of control point, side number of edge control point, 
control point index along an edge, edge control point characterization, 
global coordinates of control point, and upper surface normal at control 
point. 


6.3.4 Edge Control Point Data (OPTIONAL, IEDGEP=1) 


Edge control point data, shown in Figure 45, are printed when IEDGEP=1 
(Card 19). These data list the edge control point index, panel number of 
control point, influencing panel number, side of influencing panel, and 
global coordinates of control point. 

6.3.5 Singularity Grid Data (OPTIONAL, ISINGS=1) 

Singularity grid data, shown in Figure 46, is printed when ISINGS=1 (Card 
19). The program prints out the source strength, the doublet strength, 
and the components of the surface vorticity vector along the row and 
column directions and in direction normal thereto as well as along the 
global coordinate axis. This is done for the nine canonical points: the 
four corners, the mid points of each side, and the panel center. These 
data are useful in checking the continuity of the doublet strength from 
one network to another. Data are organized by network. 

6.3.6 Elapsed CPU Time (OPTIONAL, IPTIME=1) 
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When IPTIME=1 (Card 18), the program prints the elapsed CPU time from 
various programs and subroutines. 



GEOHCTRV DATA 


HESH POINT DATA 


NUMBER 

ROU 

COLUMN 

NET. NO. 

X 

Y 

z 

1 

1 

1 

1 

O.DOCOOCOOOO 

O.OCDOCOOOOO 

0.0000000000 

2 

j> 

1 

1 

D.OOinOOOCOD 

O.nCODODOOOO 

o.oocorocooo 

3 

3 

1 

1 

o.ooocrooooo 

O.CCOOODGOOO 

o.cooooococo 

4 

4 

1 

1 

0 .OOOOQCOCOQ 

O.ODCQOQSOOO 

0.0000000000 

5 

5 

1 

1 

O.DOCOOOOOOO 

0.0 COOOOCOCiD 

0.0003000030 

f, 

6 

1 

1 

O.ODOOOODCDO 

O.ODCODOOOOO 

0.000000 CODO 

7 

7 

1 

1 

o.DorcnoDCto 

0.0 DCC OOCOOO 

0.0000000000 

8 

1 

2 

1 

rPDCDOOCOOO 

O.OCDGOOCOGO 

0.0000000000 

9 

2 

2 

1 

.2CC3DOOOOO 

.01725000CO 

0.0000000000 

ID 

3 

2 

1 

.2CDDOCOOOC 

.0287500000 

O.OOOOOOOOGO 

11 

A 

2 

1 

.200COOODQO 

.0359375000 

O.OOOOOOCOOO 

12 

5 

2 

1 

,200Df OCCOO 

.0 A3125l'0P0 

O.OOOOOOODOD 

15 

6 

2 

1 

.20000000DO 

.0503125000 

0.0000000000 

lA 

7 

2 

1 

.2C00C0CC00 

.0575000000 

0.0000000000 

15 

1 

3 

1 

. AOODOGODCO 

o.oicooooooo 

O.DOCOOOOOOO 

16 

2 

3 

1 

.4COPDOOOCO 

.0345000000 

O.OOOODOOCOO 


FIGURE i(2 MESH POINT DATA 



SINGULARITT DISTRIBUTION DEFINITION 


PANEL DISTRIBUTION QUANTITIES 


IS - singularity parameter no. 

A1...A9 - COEFFICIENTS FOR THE 9 CANONICAL PANEL LOCATIONS 


PANEL 

M 

O 

network no. = 

1 DISTRIBUTION ORDER 

= 2 

NUMBER OF SINGULARITY PARAKETERS = 12 


IS 

A1 

A2 

A3 

A4 

AS 

A6 

A7 

AB 

A9 

1 

.10000E«01 

-.4in9fiE-Pi 

-,EIiC6''r-01 

.10000E401 

0. 

-.85571E-01 

-.2913BF-1G 

.100DOE401 

0. 

2 

0. 

.452rbE*tO 

-.b4;.4rE-02 

n. 

.1DOOOE4D1 

.4972RE-C2 

-.269D9r>02 

0. 

0. 

3 

0. 

0. 

.i4ttpr*oo 

0. 

0. 

,53159E+00 

.40235E4Pa 

0. 

.10Q00C401 

4 

c. 

n. 

.3<fP£flE40P 

0. 

0. 

-.1 1156E-02 

.60335E*U0 

0. 

0. 

b 

0. 

0. 

-.F0603E-02 

0. 

c. 

0. 

-.300G3E-C2 

c. 

0. 

10 

0. 

.69&i3E40Q 

-.4ee?2E-oi 

0. 

0. 

.27273E-01 

-.128061-01 

0. 

0. 

n 

c. 

P. 

.34b04E4t0 

0. 

0. 

.S9fc67E*00 

.12171E-C1 

0. 

0. 

12 

0. 

0. 

.4327?E400 

0. 

0. 

.19097E-01 

.13435E-ri 

0. 

0. 

13 

0. 

0. 

-.54543E-01 

0. 

0. 

Or 

-.12799E-01 

0. 

0. 

le 

D. 

-.10959E400 

0. 

0. 

0. 

-.1O017E-O1 

0. 

0. 

0. 

19 

0. 

D. 

-.51613E* .1 

0. 

0. 

-.76C63E-01 

0. 

0. 

c. 

20 

0. 

D. 

-.51H3E-01 

0. 

D. 

-.68302E-02 

0. 

0. 

0. 


PANEL 

NO. s 2 

NETWORK NO. = 

1 DISTRIBUTION ORDER 

= 2 1 

NUMBER OF SINGULARITY PARAMETERS * 14 


IS 

A1 

A2 

A3 

A4 

Ab 

A6 

A7 

A8 

A 9 

1 

»1000eE40! 

-.5H06FE-01 

-,f 8n65E-01 

. 10000E401 

-.291 38E-16 

-.65425E-01 

-.30840E-14 

.10000E401 

Cl. 

2 

0. 

b4246CT02 

0, 

0. 

-.26909E-02 

0. 

0. 

0. 

0. 

3 

0. 

. 14G6PE4G0 

-.26362C-02 

c. 

.4023SE4aO 

-.49277E-03 

-.192b5E-02 

0. 

r. 

4 

0. 

.34fi6BE400 

.132b4E400 

0. 

.60335E400 

.53091E4C0 

.3B899E400 

0. 

. 10000£«01 

5 

0. 

-.60F C3C-02 

.354P9E400 

0. 

-♦30063E-C2 

.47507E-02 

.61361E400 

0. 

c. 

6 

0. 

0. 

-.91Bt9C-C3 

0. 

0. 

0. 

-.67100E-03 

0. 

0. 

10 

c. 

-.468£2E-01 

0. 

0. 

-,126I)6E-01 

0. 

0. 

0* 

0. 

1 1 

0. 

.34P04E400 

-.2S726Er01 

0. 

.12171C-01 

.1989fiE-01 

-,4497uE-C2 

0. 

c. 

12 

0. 

.4327bE400 

.35412E400 

0. 

,13455E-01 

.5S63eE400 

.70071E-02 

0. 

1. 

13 

0. 

-.54S43E-01 

,3bS30E*00 

0. 

••.12799E-01 

.26R4CE-01 

.32207E-02 

0. 

0* 

14 

0. 

0. 

-.B 2 f,erE -02 

0. 

0. 

0. 

-.57308E-02 

0. 

0. 

19 

0. 

-.B1613E-01 

0. 

0. 

0. 

-.71435E-02 

0. 

0. 

c. 

20 

0. 

-,bl613E-Dl 

-.51613E-01 

0. 

0. 

-.75934E-01 

0. 

0. 

0. 

21 

0. 

0. 

-.5U13E-01 

0. 

0. 

-.9fi8B3E-02 

0. 

0. 

0. 


FIGURE SINGULARITY DISTRIBUTION DEFINITION DATA 



CONTROL POINT DATA 


CONTROL POINT LOCATIONS AND NORMALS 


JCN - CONTROL POINT INDEX (INCLUDING THOSE USED FOR FORCE CALCULATION 
KC - NETWORK NO. OF CONTROL POINT 
IPC - PANEL NO. OF CONTROL POINT 
ISC - SIDE NO. OF EDGE CONTROL POINT 
IZC - CONTROL POINT INDEX ALONG AN EDGE 
ICH - EDGE CONTROL POINT CHARACTERIZATION 
=0 REAL 

=1 doublet value matching 

=2 DOUBLET NORMAL DERIVATIVE HATCHING 

DOUBLET TANGENTIAL DERIVATIVE HATCHING 
X.Y.Z - GLOBAL COORDINATES OF CONTROL POINT 
NX«NV«NZ - UPPER SURFACE NORMAL AT CONTROL POINT (IN GLOBAL COORDINATES) 


JCN 

KC 

IPC 

ISC 

12C 

ICH 

X 

V 

Z 

NX 

NY 

NZ 

1 

1 

6 

3 

7 

1 

-.ooocoo 

O.CCODOC 

o.ocooco 

0.000000 

O.ODOOOD 

1 . oooooo 

2 

1 

1 

1 

? 

2 

. lOCOGO 

0 . ooocoo 

0.000000 

0.000000 

0.000000 

l.OOOQOO 

3 

1 

1 

0 

0 

0 

.lOOPDO 

.OOA3I3 

0 . ooocoo 

0.000000 

0 . oooooo 

1.000000 

A 

1 

2 

0 

0 

0 

.100000 

. OllEOO 

0.000000 

0 . oooooo 

O.COOGOO 

l.CODCOO 

5 

1 

3 

0 

0 

0 

.100000 

.018172 

0,000000 

c. 000000 

o.ooooco 

1. EOOOOO 

0 

1 

A 

0 

0 

0 

.100(100 

.019766 

0.000000 

0.300000 

0.000000 

1. OOOOOO 

7 

1 

5 

0 

G 

0 

.100000 

. 023355 

o.ooooco 

0.000000 

0.000000 

1.000000 

B 

1 

G 

p 

0 

0 

.100000 

.026953 

0 . ooocoo 

0.000000 

0 . oooooo 

l.CDDOOO 

a 

1 

6 

3 

b 

2 

.looroo 

.02f 750 

0 . ooocoo 

0.000000 

0.000000 

l.CCCOOO 

K 

1 

7 

1 

3 

2 

.30C0C0 

O.OCnCCiO 

o.ooooco 

O.OOO'JOO 

o.ooooco 

1 . oooooo 

11 

1 

7 

0 

0 

0 

.30CC00 

.012538 

o.oocroo 

o.ocoooo 

0.300000 

i.tooooo 

12 

1 

8 

0 

0 

0 

.300000 

.03A500 

o.oocooo 

0.000000 

0 . oooooo 

1.000000 

13 

1 

9 

0 

0 

0 

.300000 

. 0A851 6 

0.000000 

0.000000 

o.ooooco 

1.000000 

J A 

1 

1C 

0 

0 

0 

.300000 

.059297 

0.000000 

0.000000 

0 . oooooo 

l.CQbOGO 

15 

1 

11 

0 

0 

0 

.300000 

.070078 

O.OCOOOD 

0.000000 

o.ocooco 

l.FODOOO 

16 

1 

12 

0 

0 

Q 

.300000 

.080659 

o.ooocoo 

c. oooooo 

0 . oooooo 

1 . oooooo 

17 

1 

12 

3 

5 

2 

.30CCOO 

.086250 

0.000000 

0.000000 

0.000000 

i.cocooo 


192 

7 

117 

0 

0 

0 

1.O75P0O 

.212122 

.232052 

.OOOAOO 

.970317 

-.2A1P35 

193 

7 

110 

0 

0 

0 

1.075000 

.198191 

.175756 

-.010926 

.970260 

-.2A1B20 

19A 

7 

IIB 

3 

2 

0 

1.075000 

.191075 

.1A7607 

-.016586 

.97016A 

-.2A1802 


NF (NUMBER OF SINGULARITY STRENGTH PARAMETERS) = 173 

NO (NUMBER OF PANEL ORIENTATION ANCLES) = b6 

QQ NH (NUMBER OF GEOMETRY PARAMETERS - LAMBDA AND NUI = lA 

VO 


FIGURE CONTROL POINT DATA 



tDGt CONTROL POINT OfcTik 


JC - EDGE CONTROL POINT INDEX 
IPC - PANEL NO. OF CONTROL POINT 
IP - INFLUENCING PANEL NO. 

IS - SIDE NO. OF INFLUENCING PANEL 
ZX,ZT.ZZ - GLOBAL COORDINATES OF CONTROL POINT 


JC 

IPC 

IP 

IS 

zx 

ZY 


2Z 

2 

1 

1 

1 

•100C0OOOE«OO 

0. 

0. 


JC 

IPC 

IP 

IS 

zx 

ZY 


Z2 

1 

6 

6 

3 

-.99908921E-15 

0. 

0. 


JC 

IPC 

IP 

IS 

ZX 

ZY 


ZZ 

9 

6 

6 

3 

.lOOOOODCE^On 

.267500DDE-01 

0. 


JC 

IPC 

IP 

IS 

ZX 

ZY 


22 

10 

7 

7 

1 

.300nODOOE-fOO 

0. 

Q. 


JC 

IPC 

IP 

IS 

zx 

ZY 


ZZ 

1 

6 

A9 

1 

-.994CR921E-15 

0. 

0. 


JC 

IPC 

IP 

IS 

ZX 

ZY 


ZZ 

9 

6 

99 

1 

.lRODCOCOE+00 

.2B750000E-01 

0. 


JC 

IPC 

IP 

IS 

ZX 

ZY 


ZZ 

17 

12 

12 

3 

.300000POE*03 

.e6250000E-Cl 

0. 


JC 

IPC 

IP 

IS 

ZX 

ZY 


ZZ 

16 

13 

13 

1 

.500II00(I0E*C0 

0. 

0. 


JC 

IPC 

IP 

IS 

ZX 

ZY 


ZZ 

25 

la 

le 

3 

.500000COE+aO 

.143750COE*00 

0. 


JC 

IPC 

IP 

IS 

zx 

ZY 


ZZ 

26 

19 

19 

1 

»7ooooocor«oo 

0. 

0. 


JC 

IPC 

IP 

IS 

2K 

ZY 


ZZ 

17 

12 

57 

1 

.300000CCE>00 

.B6250Q00E>C1 

0. 



FIGURE ‘\5 EDGE CONTROL POINT DATA 



IP - P*NtL NO. 

1»J - ROW ANll COLUMN INDICPS OF THF 9 CANONICAL POINTS 
- GLOBAL COCPDlNAirS OF THE flVEN POINT 
SO - SOURCE SINGULARITY STRENGTH 
00 - nOUBLET SINGULARITY STRENGTH 
DAfOYtOZ - VORTICITY COMPONENTS IN X,Y,Z DIRECTIONS 
OMfON - VCiRTICITY COMPONENTS IN M,N DIRECTIONS 
DMP.DNP - VOPTICITY COUPONENTS IN DIRECTIONS PERPENDICULAR TO H«N 


NETWORK NO. = I NETWORK TYPE - 2 


SINGULARITY GRID 


1? 

I 

J 

X 

Y 

Z 

so 

DO 

DX 

DY 

D2 

DM 

DN 

PHP 

DNP 

I 

1 

1 

G.bOODO 

0.00000 

O.OGODC 

C. 00000 

•00000 

-.cnooo 

.30504 

0. 00000 

0.00000 

-.00000 

0.00000 

.30504 

1 

2 

1 

O.OOCOO 

0.00000 

0. CbOOD 

0.00000 

.00500 

-.OOCOO 

.30340 

0.50000 

0.00000 

.01307 

0. OCOOO 

•30312 

1 

3 

1 

0. DOOl'O 

o.coooo 

P. OOCOO 

O.DOODO 

.ocooo 

- .oorco 

.29602 

O.OODFO 

O.OOOOP 

.02561 

0.00000 

•29692 

1 

1 

2 

. ICOPO 

o.ropoo 

O.COOOO 

0.00000 

.0300 A 

O.OODOD 

.29579 

O.OOOOD 

.29579 

r.PODOO 

-.00000 

.25579 

1 

2 

2 

.10000 

.00<)31 

0.00000 

0.00000 

.029RB 

.07356 

. 29744 

0.00000 

.29744 

.08632 

-.07358 

.29399 

1 

3 

2 

.iCiono 

.bOR63 

C..OCOUO 

0.00000 

.02^41 

.14715 

.30281 

0.00000 

.30281 

.17263 

-.14715 

.23905 

1 

1 

3 

.20000 

0.00000 

\). 00000 

0.00000 

. 05916 

.00507 

• 2B653 

0.00000 

.28653 

•00507 

-.00507 

•28653 

1 

2 

3 

.2'>000 

.008F3 

0.00000 

o.onooo 

.05865 

.06574 

.28796 

0.00000 

.28796 

.07609 

-.06574 

.28466 

1 

3 

3 

.20000 

.01 725 

0 . ocooo 

O.OCCGO 

.058C2 

.12642 

.29312 

D.GOOCO 

.29312 

.151 14 

-.12642 

.28117 


2 

1 

1 

0. OCOCO 

-.00000 

0.00000 

0.05000 

. 00500 

.OOCOO 

.29802 

o.coooo 

0.00000 

.02561 

0.00000 

.29692 

2 

2 

1 

O.COOOO 

o.coooo 

0. OCOOO 

O.POOOO 

. 00000 

.OOCOO 

.29457 

0.00000 

O.OOCOO 

• 03365 

O.OOOOD 

.29264 

2 

3 

1 

0. OCOOO 

-.00000 

0.00000 

*0.00000 

.00000 

.CODOO 

.26807 

0.00000 

O.OOCOO 

.04099 

0.00000 

.28514 

2 

1 

2 

.10550 

.C0662 

O.OOCOO 

0.00000 

.02941 

.08804 

.29771 

0.00000 

.29771 

.11330 

-.08804 

.28905 

2 

2 

2 

.19090 

.01150 

O.DCPOD 

o.oonoo 

.02904 

.16996 

.30575 

0.00000 

,30570 

.20376 

-.16996 

•28428 

c 

3 

2 

.10000 

.01 438 

D. OCOOO 

0.00000 

.02643 

.25189 

.31673 

O.DOODO 

.31673 

.29439 

-.25189 

•277C7 

2 

1 

3 

.20000 

.01 725 

P.COOUO 

0. DOOOO 

.05802 

.10867 

.29159 

O.OOOOD 

.29159 

.13332 

-.10867 

.28117 

2 

2 

3 

.20000 

.02300 

O.COOOO 

0.00000 

. 05723 

.16686 

.29694 

O.COOOO 

.29694 

.19969 

-.16686 

•27593 

2 

3 

3 

.20000 

.02875 

O.OOCOO 

0.00000 

. 05610 

.22505 

.30533 

0.00000 

.30533 

.26620 

-.22505 

•27020 


3 

1 

1 

-.05000 

-.00050 

0. OOFOO 

C. 00000 

.OOCOO 

-.OOCOO 

.28807 

0.00000 

0 . ocooo 

.04059 

O.ODOCO 

•28514 

3 

2 

1 

0. OCOOO 

0.00000 

Q.DCPOO 

O.OOCOO 

• ocooo 

-.OOCOO 

•28408 

0.00005 

C. 00000 

.04535 

0 . ocooo 

.28044 

3 

3 

1 

-.OOCOO 

-.OOCOO 

n. OCOOO 

O.OOCOO 

.00000 

-.OODCO 

.27794 

O.OOCOO 

0.00000 

.04915 

0 . ocooo 

.27355 

3 

1 

2 

.10000 

.01438 

o.coooo 

n.pccoc 

.02843 

.18064 

.30649 

0.00000 

.30649 

.22241 

-.18C64 

•277E7 

3 

2 

2 

. JOOOO 

.01617 

0. OODDO 

O.OOCOO 

• 0280? 

.27310 

.32053 

0.00000 

.32053 

.32077 

-.2731C 

•27281 

3 

3 

2 

.10000 

.01797 

0.00000 

O.OOQOC 

.02745 

.36556 

.33672 

o.ooooc 

.33672 

.41935 

-.36556 

•26676 

3 

1 

3 

.25005 

.C2R75 

O.COOOO 

0.C50C0 

.05615 

.19940 

.30164 

0.00000 

.30164 

.24C29 

-.19940 

.27020 

3 

2 

3 

.20050 

.03234 

0. OIJOOD 

0.05000 

.05527 

.26417 

.31136 

0.00000 

.31136 
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FIGURE 46 SINGULARITY GRID DATA 


6.3.7 Near Field/Far Field Information (OPTIONAL, IPNPIC=1) 


When IPNPIC=1 (Card 18) program VINFCC prints a table of panel influence 
coefficients count for source and doublet. The counts are listed 
separately for no influence, monopole far field, dipole far field, 
quadrupole far field, one subpanel intermediate field, two subpanel 
intermediate field, and eight subpanel near field. 

6.3.8 Out-of-Core Solver Information (OPTIONAL, IPS0LV=1) 

When IPS0LV=1 (Card 18), the PTSOLV solution package prints a description 
of the problem being solved and an error analysis report. 

6.4 TAPE14 Save File 


A data file is set up on TAPE14 which can be used for post processing or 
for a solution restart. Network mesh point data is generated during a 
data check which can be used for post processing graphics. The user must 
provide his own interface and graphics software. During the iterative 
solution network and singularity data necessary for a solution restart are 
saved. Section 7.3 describes the file usage. The TAPE14 format is given 
in Table 5. 


92 



TABLE 4 DEFINITION OF OUTPUT QUANTITIES 


QUANTITY 


DEFINITION 


WHEN 

PRINTED 


JC 

IP 

(X,Y,Z) 

DO 

(DX,DY,DZ) 

SO 

(FSVX,FSVY, 

FSVZ) 

PHIU 

(WXU,WYU, 

WZU) 

PHEU 

(PWXU,PWYU, 

PWZU) 

CPLINU 

CPSLNU 

CPSLNDU 

CPI2NU 

PHIL 

(WXL.WYL, 

WZL) 

PHEL 

(PWXL,PWYL, 

PWZL) 

CPLINL 

CPSLNL 

CP2NDL 

CPISNL 

WNU 

WNL 

WTU 

WTL 

PHIUI 

PHILI 

PWNU 

PWNL 


Physical Quantities 
Cumulative control point index 
Index of panel containing control point 
Global coordinates of control point 
Doublet strength 

Global Coordinates of surface vorticity 
vector 

Source strength 

Freestream velocity vector in global 
coordi nates 

Upper surface total potential 
Upper surface total mass flux vector in 
global coordinates 

Upper surface perturbation potential 
Upper surface perturbation mass flux vector 
in global coordinates 
Upper surface linearized pressure coefficient 
Upper surface slender body pressure coefficient 
Upper surface second order pressure coefficient 
Upper surface isentropic pressure coefficient 
Lower surface total potential 
Lower surface totol mass flux vector in global 
coordinates 

Lower surface perturbation potential 
Lower surface perturbation mass flux vector 
in global coordinates 

Lower surface linearized pressure coefficient 
Lower surface slender body pressure coefficient 
Lower surface second order pressure coefficient 
Lower surface isentropic pressure coefficient 
Normal component of upper surface total mass 
flux vector 

Normal component of lower surface total mass 
flux vector 

Magnitude of tangential component of upper sur- 
face total mass flux vector 
Magnitude of tangential component of lower sur- 
face total mass flux vector 
Upper surface total mass flux potential 
Lower surface total mass flux potential 
Normal component of upper surface perturbation 
mass flux vector 

Normal component of lower surface perturbation 
mass flux vector 


Every ITPRINT 
iteration, 
initial 
solution, 
final iteration 
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TABLE 4 (CONTINUED) 


QUANTITY 

CPLIND 

CPSLND 

CP2NDO 

CPISND 


AREA 

(FX,FY,FZ) 

(MX,MU,MZ) 


NET 

NUP 

NATl 

NAT2 

NSDl 

NSD2 

NCRl 

lEDGAl 

IEDGA2 

KEDGl 

KEOG2 

MEDGIT 

MEDGIF 

MEDG2I 


DEFINITION 


WHEN 

PRINTED 


Difference between upper and lower surface 
linearized pressure coefficient 
Difference between upper and lower surface 
slender body pressure coefficient 
Difference between upper and lower surface 
second order pressure coefficient 
Difference between upper and lower surface 
isentropic pressure coefficient 


Forces and Moments Data 


Always 


Total area of panels 
Global coordinates of force coefficient 
(Upper surface, lower surface, difference) 
Moment coefficients about global principal 
axes (Upper surface, lower surface, differ- 
ence) 


Update Index Arrays 


Always 


Network sequence number 
Network update index 

Square number of network to which side 1 
of NET is attached 

Sequence number of network to which side 4 
of NET is attached 

Side of NATl to which side 1 of NET is 
attached 

Side of NAT2 to which side 4 of NET is 
attached 

Leading corner of NSDl to which side 1 of NET 
is attached 

Cumulative index of points on side 1 which 
have been assigned matching points 
Cumulative index of points on side 4 which 
have been assigned matching points 
Same definition as NATl when NATl = 0 
Same definition as NAT2 when NAT2 = 0 
Row index of point on NATl to which initial 
point on side 1 of NET is attached 
Row index of point on NATl to which final 
point on side 1 of NET is attached 
Row index of point on NAT2 to which initial 
point on side 4 of NET is attached 
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TABLE 4 (CONCLUDED) 


QUANTITY 


DEFINITION 


WHEN 

PRINTED 


MEDG2F 

MED61I 

NEDGIF 

NEOG2I 

NED62F 


Row index of point on NAT2 to which final 
point on side 4 of NET is attached 
Column index of point on NATl to which 
initial point on side 1 of NET is attached 
Column index of point on NATl to which final 
point and side 1 of NET is attached 
Column index of point on NAT2 to which initial 
point and side 4 of NET is attached 
Column index of point on NAT2 to which final 
point and side 4 of NET is attached 


S 

IP 

I 

J 

(X,Y,Z) 

SO 

DO 

(DX,DY,DZ) 

(SM,SN) 

(SMP,SNP) 


Singularity Grid Data (ISINGS = 1) 

Solution singularity parameters 
Cumulative index of panel on which 
singularity distribution is evaluated 
Local row index of evaluation point 
Local column index of evaluation point 
Global coordinates of evaluation point 
Source strength value at evaluation point 
Doublet strength value at evaluation point 
Global coordinates of surface vorticity 
vector at evaluation point 
Derivative of doublet strength in (row, 
column) directions respectively 
Derivative of doublet strength in directions 
normal to (row, column) directions respect- 
ively 


TABLE 5 TAPE14 FORMAT 


(1) When 3DATA CHECK is specified, the following FORTRAN statements in program 
INPUT are used to save network mesh points on NSAV=TAPE14 (format 6F10.5): 

FNT = NNEH 

WRITE (NSAV, 1030) FNT 
1030 FORMAT (6F10.5) 

DO 1070 K = 1, NNETT 
NMK = NM(K) 

NNK = NN(K) 

FM = NMK 
FN = NNK 

WRITE (NSAV, 1030) FM,FN 

DO 1060 J = 1, NNK 

JM = (J-1) * NMK + NZA(K) 

WRITE (NSAV, 1030) (ZM(1,I+JM), ZM(2,I+JM), ZM(3,I+JM), I = 1, NMK) 
1060 CONTINUE 
1070 CONTINUE 

(2) When ^ITERATION or ^LEAST SQUARES ITERATION is specified, the network mesh 
points and values of singularity parameters for the current iteration (no. 
JT) are saved on NSAV=TAPE14 (unformatted) using the following FORTRAN 
statements : 

NNETPl = NNETT + 1 
REWIND NSAV 
WRITE (NSAV) JT 

WRITE (NSAV) NNETT, NSNGT, NZMPT, NNETPl 
WRITE (NSAV) (NZA(I), I = 1, NNETPl) 

WRITE (NSAV) (ZM(1,J), ZM(2,J), ZM(3,J), J = 1, NZMPT) 

WRITE (NSAV) (S(I), I = 1, NSNGT) 
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7.0 COMPUTER PROGRAM DESCRIPTION 


This computer program is written in the CDC FORTRAN Extended (FTN4) language 
for the CDC Network Operating System (NOS). It uses overlay structures and 
fourteen disk files which include the standard system files INPUT (TAPES) for 
card reading and OUTPUT (TAPES) for printing. The program has been checked 
out and run on the Langley Research Center's CDC CYBER series computers. 

The computer code implements recent advances in the solution of 
three-dimensional flow over wings with leading edge vortex separation. It has 
been designed and developed for the purpose of performing numerical experiment 
studies with the flow model. 

The code includes two iterative solution procedures: (i) Quasi-Newton scheme 

and (ii) Least Squares Method. The least squares procedure for damping vortex 
sheet geometry update instabilities was developed to alleviate the convergence 
problem for certain cases using the Quasi-Newton ' iterative scheme. It is 
restricted to run smaller problems (see discussion in User's Input Guide) in 
the present set up and takes more computational time to execute. In the 
future, we hope to further develop the least squares procedure so that it can 
be used to execute larger problems as well as taking less computational time. 

7.1 Basic Program Structure 

The computer program consists of one main overlay, six primary overlays, 
three secondary overlays and one user library. A schematic diagram of 
basic program structure is illustrated in Figure 47. 

7.2 Description of Overlay Programs 

The following is a discussion of the overlay programs. A detailed flow 
chart of the main overlay program A378 is shown in Figure 48. 

7.2.1 OVERLAY (MAIN, 0,0) 

Main OVERLAY (MAIN, 0,0) 

Program A378 

Purpose To perform various tasks by calling the following overlay 
programs and subroutines: 

0 Program INPUT to process the input data and set up network 
mesh points 

0 Program CONFIG to compute panel geometry, panel singularity 
distribution and panel control points defining quantities 

0 Program VINFCC to calculate and store induced potential and 
velocity coefficients 
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FIGURE i)7 BASIC PROGRAM STRUCTURE 













FIGURE 48 FLOW CHART OF MAIN OVERLAY PROGRAM A378 
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FIGURE 48 CONTINUED 


100 









FIGURE 48 CONCLUDED 
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Discussion 


Subroutine INIMU to compute initial values for singularity 
strength parameters 

Program OUTPUT to calculate and print the singularity grid 
and the physical flow quantities for each network 

Subroutine JTFLOW to perform subiteration for a solution of 
singularity strength parameters 

Subroutine ITFLOW using Quasi-Newton scheme to find an 
iterative solution to the flow problem 

Subroutine LSFLOW using Least Squares Method to find an 
iterative solution to the flow problem 

The main overlay program A378 first sets up data blocks for 
transferring among the overlay programs and also 
initializes the data. At the beginning of the execution of 
the code, disk files except system files INPUT (TAPES), 
OUTPUT (TAPE6), and random access file (TAPE4) are 
rewound. The code then reads and prints a list of the 
input data cards. Program INPUT is first called to process 
the input data and set up network mesh points. Next 
program CONFIG is called to compute panel geometry, panel 
singularity distribution and panel control points defining 
quantities by using the secondary overlay programs TGEOMC, 
TSING and TCNTRL. Induced potential and velocity influence 
coefficients, are then calculated via program VINFCC. 

If the input data indicates no previous iteration results 
are to be used, then the code calls subroutine INIMU to 
calculate the initial singularity strength parameters 
values. Otherwise, the code reads in the values provided 
by the user on disk file TAPE14. If no iteration is 
requested, then program OUTPUT is called to calculate and 
print the singularity grid, physical flow quantities, and 
mesh points for each network. 

If the user requests only the subiteration then the code 
calls subroutine JTFLOW for an iterative solution of the 
singularity strength parameters values, and program OUTPUT 
for printing the results of physical flow quantities and 
network mesh points. 

When full iteration using Quasi -Newton scheme or least 
squares method is requested, then the code will proceed 
after the subiteration being exercised to find an iterative 
solution to the flow problem by using either subroutine 
ITFLOW or subroutine LSFLOW. 
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7.2.2 OVERLAY (DATA, 1,0) 


Primary 

Program 

Purpose 


Discuss 


OVERLAY (DATA, 1,0) 

INPUT 

To read and process the input data as follows: 

0 Set up network mesh points by using various preprocessors 

0 Calculate free stream velocity and compressibility 
direction and metric A and B 

o Determine all network edge abutments and abutment 
intersections 

0 Obtain initial panel orientation angles and geometry 
parameters X and v 

ion The input data cards are processed as indicated in the 

discussion of the User's Input Guide. After reading in the 
physical quantities such as angle of attack, yaw angle, 
symmetry or asymmetry, Mach number, and reference values, 
the code checks which of the following option is 
requested: (1) data check, (2) subiteration, (3) full 

iteration using Quasi-Newton scheme, or (4) full iteration 
using Least Squares Method. The printing options are read 
in next. Finally, the code reads the number of networks 
and the specified preprocessor. The preprocessors such as 
^QUADRILATERAL, ^GOTHIC, 2V0RTEX, and ^TRAILING WAKE are 
then called to set up the network mesh points. The code 
proceeds to calculate free stream velocity, compressibility 
direction matrices A and B, and orthogonal transformation 
matrix from reference coordinates into wing axis 
coordinates. If previous iteration results are to be used, 
network mesh points data will be read in from disk file 
TAPE14. On the other hand, if 3DATA CHECK is specified, 
network mesh points data will be saved on disk file NSAV 
for external graphic processing. A printer plot of cross 
sections of the initial vortex system will be produced when 
this option is chosen. A printout of input network mesh 
points data is always provided by the program. Before 
returning the control to the main overlay, the code calls 
subroutine ABTCAL to determine all network edge abutments, 
and abutment intersections, and also obtains initial panel 
orientation angles and geometry parameters and 
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7.2.3 OVERLAY (NETGCS, 2,0) 


Primary OVERLAY (NETGCS. 2,0) 

Program CONFIG 

Purpose To compute panel geometry, panel singularity distribution, 

and panel control points defining quantities for all 

networks . 

Discussion This overlay program serves as a driver for calling the 
following three secondary overlay programs: 

0 Program TGEOMC to compute panel geometry defining 
quantities 

0 Program TSING to compute panel singularity 
distribution defining quantities 

0 Program TCNTRL to compute panel control points 
defining quantities 

7.2.4 OVERLAY (NETGCS, 2,1) 

Secondary OVERLAY (NETGCS, 2,1) 

Program TGEOMC 

Purpose To compute panel geometry defining quantities for all 

networks 

Discussion If diagnostic geometry information is desired (IGEOMP = 1), 
the code prints mesh points data along with their row and 
column indices of the network. The main function of the 
code is to call subroutine GEOMC for each network to 
compute panel geometry defining quantities. It also sets 
up arrays containing number of panels in each network 
(NP(K)), number of mesh points in each network (NZ(K)) and 
accumulated sum of NP (NPA(K)), and accumulated sum of NZ 
(NZA(K)); and obtains the total number of panels (NPANT) 
and the total number of mesh points (NZMPT) for all 
networks . 
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7.2.5 OVERLAY (NET6CS, 2 , 2 ) 

Secondary OVERLAY (NETGCS, 2.2) 

Program TSING 

Purpose To compute panel singularity distribution defining 

quantities for all networks 

Discussion The code calls various routines depending on the type of 
each network to compute panel singularity distribution 
defining quantities: 

0 Routine SING for type 1 network 

0 Routine DASPL for type 2 network 

o Routine DDSPL for type 4 network 

0 Routine DSDSPL for type 6 network 

0 Routine DWSPL for type 8, 10, 14, 16 network 

It also sets up arrays containing number of singularity 
parameters in each network (NS(K)), and accumulated sum of 
NS(NSA(K)), and obtains the total number of singularity 
parameters (NSNGT) for all networks. 

7.2.6 OVERLAY (NETGCS, 2,3) 

Secondary OVERLAY (NETGCS, 2,3) 

Program TCNTRL 

Purpose To compute panel control points defining quantities for all 

networks 

Discussion The code calls subroutine CONTRL for each network to 

compute panel control points defining quantities. It also 
sets up arrays containing number of control points in each 
network (NC(K)) and accumulated sum of NC (NCS(K)), and 
obtains the total number of control points (NCTRT) for all 
networks. Next, defining quantities of some special 
control points used for calculating fed sheet force are 
computed by calling subroutine CONFRC. The sum of the 
total number of these special control points and that of 
the original control points is given as a parameter NCTRTE. 
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7.2.7 OVERLAY (AICGEN, 3,0) 


Primary 

Program 

Purpose 

Discussion 


OVERLAY (AICGEN, 3,0) 

VINFCC 

To calculate and store induced potential and velocity 
coefficients 

The code first obtains far field moments for hyperboloidal 
panels by calling subroutine FFHPM6. It starts to compute 
the potential /velocity influence coefficients by storing 
defining quantities for a group of control points in the 
available core of a scratch array. Then it proceeds to 
calculate panel influence coefficients by looping through 
all panels for that group of control points. After these 
calculated influence coefficients are being stored in a 
random access disk file, the code goes back to store 
defining quantities for another group of control points and 
perform the same calculation as described before. This 
process terminates when it is done with all control 
points. At the very first call of this overlay, all of the 
influence coefficients wil be computed and stored in a 
random access disk file. After that, only the part of the 
influence coefficients affected by the perturbation of 
geometry will be calculated and replaced on the random 
access disk file. This cost saving scheme is controlled by 
a parameter NRAIC passing through the common block REAIC. 

The potential /velocity influence coefficients of those 
special control points used for fed sheet force calculation 
are also obtained in this overlay program via subroutine 
AICFOR. Finally, information for far field, intermediate 
field, and near field are provided when requested by the 
user (IPNPIC=1). 



7,2.8 OVERLAY (JACGEN, 4,0) 


Primary 
Pro gr am 
Purpose 
Discussion 


OVERLAY (JACGEN, 4,0) 

AJGEN 

To calculate and store the analytic Jacobian matrix. 

The code first sets the row and column dimensions of the 
Jacobian matrix according to the type of iteration 
requested: 

(i) Subiteration (NDZA=0) - iteration performed only on 
the singularity strength parameters 

(ii) Iteration without force boundary condition (NDZA=1) - 
iteration performed on the singularity strength 
parameters and panel orientation angles 

(iii) Full iteration (N0ZA=2) - iteration performed on the 
singularity strength parameters, panel orientation 
angles, and geometry parameters A and v. 

(iv) Least Squares iteration (NDZA=3) - iteration performed 
on all the parameters with additional equations 
consisting of twist function 

If it is not subiteration (NDZA=0), a subroutine ZTHET is 
called to calculate and store perturbation of network mesh 
points with respect to panel orientation angle and geometry 
parameters A and v . The major part of the code is the 
loop that ranges over all regular control points to compute 
the Jacobian corresponding to the equations of the function 
F (see section 4.3) and store it by row on the disk file 
NSC3. For those interior control points at the network 
with update index greater than 3, it also calculates rows 
of Jacobian corresponding to the function G and saves on a 
scratch disk file if NDZA is not zero. Later, these rows 
of Jacobian are transferred and stored behind the Jacobian 
matrix corresponding to the function F. If full iteration 
is requested, the code computes rows of Jacobian 
corresponding to the function H (force boundary condition) 
via subroutine DFRC and stores also on the disk file NSC3. 

The code will compute and store rows of Jacobian 
correspondi ng to the function K (twist function) via 
subroutine DKCAL when it is desired to have least squares 
iteration. 



7.2.9 OVERLAY (SOLVER, 5,0) 


Primary 

Program 

Purpose 

Discussion 


OVERLAY (SOLVER, 5,0) 

TINVER 

To provide an interface for using a large out-of-core 
equations solver 

The input coefficient and right-hand side matrices are 
read in by row, and written into square blocks on a 
random access disk file via subroutine BLOCKR. The 
code sets up the argument list and calls the 
out-of-core equations solver PTSOLV. The solution 
matrix obtained from the solver is then unblocked and 
written by row on the disk file NANS via subroutine 
RBLOCK. 
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7.2.10 OVERLAY (RESULT, 6,0) 


Primary 

Program 

Purpose 


Discussion 


OVERLAY (RESULT, 6,0) 

OUTPUT 

To compute and print or save on a disk file the following 
quantities: 

0 Network mesh points and values of singularity strength 
parameters for each iteration 

0 Singularity grid (singularity strength and derivatives 
at 9 canonical points on each panel) on each network 
for diagnostic purpose 

0 Physical flow quantities of interest such as average, 
upper and lower surface potentials and velocities, 
singularity strength and gradient, and upper, lower 
and difference pressure coefficients, force and moment 
coefficients for each network 

0 Printer plot for cross sections of vortex systems 

At the beginning of the code, current iteration number, 
network mesh points and values of singularity strength 
parameters are saved on a disk file NSAV. These results 
could be used to restart another run later for more 
iterations. If diagnostic option (ISINGS=1) is chosen, 
singularity strength and derivatives at 9 canonical points 
on each panel are computed and printed along with the 
global coordinates of the representative points for each 
network . 

The next major portion of the code is to compute and print 
various physical flow quantities for all panel center 
control points on each network. Physical flow quantites of 
interest include average, upper and lower surface 
potentials and velocities, singularity strength and 
gradient, and upper, lower and difference pressure 
coefficients. The pressure coefficients are also stored 
for use in computing network force and moment coefficients 
by calling the subroutine FMCAL. 

Finally, a printer plot for cross sections of vortex system 
(when IPL0TP=1) and a list of network mesh points are given 
as part of the output for each iteration. 
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7,3 File Usage 


There are fourteen disk files used in the computer program. They all have 
symbolic names except TAPE4 which is used as a random access file. The 
following table shows the common block through which the disk file is 
passed, program or subroutine that uses it, and how it is being used. 


Disk 

Symbolic 

Common 

Program 

Usage 

No. 

Name 

Block 

Subroutine 


1 

NTD 

DRW I 

IDTRNS 

Store singularity spline derivatives 




OTRNS 

Retrieve singularity spline 


NTP 

PRWI 

IPTRNS 

Store panel geometry defining quantities 




PTRNS 

Retrieve panel geometry defining quantities 

2 

NTS 

SRWI 

ISTRNS 

Store singularity spline defining 
quantities 




STRNS 

Retrieve singularity spline defining 
quantities 

3 

NTC 

CRWI 

ICTRNS 

Store control points defining quantities 




CTRNS 

Retrieve control points defining quantities 

4 




Random access file declared in A378 




VINFCC, 

Store/retrieve potential and velocity 




AICFOR 

influence coefficients 




FGCAL, 

Retrieve potential and velocity 




HCAL, 

INIMU, 

AJGEN, 

etc. 

influence coefficients 

5 

NTS IN 

CM03 

A378, 

Standard system file INPUT for card 




INPUT, 

reading 

6 

NTSOUT 

CM03 

A378, 

INPUT, 

OUTPUT, 

etc. 

Standard system file OUTPUT for printing 

7 

NAIC 

SOLNT 

LSFLOW 

Store/retrieve information of 


decomposition for least squares' Jacobian 
matrix 


110 



Di sk 

Symbolic 

Common 

Program 

Usage 

No. 

Name 

Block 

Subrouti ne 


8 

NRHS 

SOLNT 

INIMU, 

Store the right-hand side matrix by row; 




ITFLOW, 

also serve as a scratch file in LSFLOW 




JTFLOW 

LSFLOW 

for temporary storage 




BLOCKR 

Retrieve the right-hand side matrix stored 
by rows for setting up square blocks to be 





used by the out-of-core equations solver 

9 

NANS 

SOLNT 

RBLOCK 

Store the solution matrix by row 




A378, 

Retrieve the solution matrix stored by 




ITFLOW, 

JTFLOW 

row 




CTRS 

Retrieve the solution matrix stored by 
row, and store it by column 




LSFLOW 

Retrieve the solution matrix stored by 
columns 

10 

NSCl 

SOLNT 


Random access file declared in TINVER 




BLOCKR 

Store square blocks of the coefficient and 
the right-hand side matrices for the 
out-of-core equations solver 




PTDCOM, 

Retrieve square blocks of the coeff- 




PTFSUB, 

cient and the right-hand side matrices. 




etc . 

and store square blocks of the solution 
matrix 




RBLOCK 

Retrieve square blocks of the solution 
matrix from the out-of-core equations 
solver 




CTRS 

Serve as a scratch file for storing matrix 
by blocks 

11 

NSC2 

SOLNT 

BLOCKR 

Random access file declared in BLOCKR 
Serve as a scratch file for setting up 
square blocks of input matrixes which are 
originally stored by row 

12 

NSC3 

SOLNT 

INIMU 

Store the coefficient matrix by row 




AJGEN 

Store the Jacobian matrix by row 




ITFLOW 

Retrieve the Jacobian matrix stored by 
row; also store the updated Jacobian 
matrix by row 




BLOCKR 

Retrieve the coefficient or Jacobian 


matrix stored by row for setting up square 
blocks to be used by the out-of-core 
equations solver 
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Disk 

Symbolic 

Common 

Program 

Usage 

No. 

Name 

Block 

Subroutine 


13 

NSCR 

SOLNT 

ITFLOW 

Serve as a scratch file for storing the 
Jacobian matrix 




LSFLOW 

Store the part of the Jacobian matrix 
corresponding to the function F 




AJGEN 

Serve as a scratch file for storing rows 
of Jacobian corresponding to the function 

G. 

14 

NSAV 

CMOS 

INPUT 

Save network mesh points when 2DATA CHECK 
is specified 




OUTPUT 

Save iteration no., part of network 
indices, network mesh points, values of 
singularity parameters 




INPUT 

Retrieve iteration no., part of network 
indices, network mesh points from previous 
iteration results 




A378 

Retrieve values of singularity parameters 


from previous iteration results 

The computer program uses the following CDC system utilities for random access 
files: 

OPENMS - declare a disk file to be random 
CLOSMS - close a random access file 
READMS - read a random access record 
WRITMS - write a random access record 
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7.4 Common Block Definition 


Variables of the more essential common blocks shared by overlay programs 
and subroutines are defined below. 


Common 

Block 

Variables 

Description 

ACASE 

ALPHA 

Angle of attack in degrees 


BETA 

Yaw angle in degrees 


FSVM 

Magnitude of free stream velocity 


FSV 

Free stream velocity vector 


OMEG 

Roll rate vector 


RC 

Center of rotation 

ADR 

RTD 

Degrees in unit radian 


DTR 

Radians in unit degree 

AICQ 

Q 

Q(I,K) = dependence of Ith Taylor's series 


B 

coefficient of panel doublet distribution 
on Kth canonical panel doublet value 
Dependence of coefficient of linear 


C 

distribution on triangle on values at 
vertices 

Dependence of coefficients of quadratic 


D 

distribution on triangle on values at 
vertices and edge lambdas 
Dependence of coefficients of reduced 


AR 

cubic distribution on triangle on values 
at vertices and center and edge lambdas 
Transformation matrix from global to local 


ARI 

coordinates 
Inverse of AR 


ARP 

Matrix transforming combined 


X 

potential /velocity vector from local to 

global coordinates 

Local coordinates of control point 


DU 

Dependence of potential /velocity vector in 


DV 

local coordinates on source distribution 
coefficients 

Dependence of potential /velocity vector in 


DVP 

local coordinates on doublet distribution 
coefficients 

Dependence of potential /velocity vector in 


PC 

global coordinates on doublet distribution 
coefficients 

Hyperboloidal panel geometry coefficient 


EN 

vectors 

Unit normal (in global coordinates) to 


P 

plane panel 

Local coordinates of panel corner points 
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Common 

Block 


ANGLEN 


BCARY 


Variables 


Description 


ZET 

ZSTP 

ZP 

R 


DS 


ZA 

ZL 

ZAP 

ZLF 

ALAM 

ANU 

NZAT 

NZAFT 


Hyperboloidal panel geometry coefficient 
vectors in local coordinates 
Parameter values corresponding to ZP 
Average of triangle corner points 
Dependence of coefficients of reduced 
cubic doublet distribution on 9 canonical 
panel doublet parameters 
Dependence of doublet strength and 
vorticity at a point on 9 canonical 
doublet parameters 

Panel orientation angles for free sheet 
network 

Panel edge length along spanwise cut for 
free sheet network 

Panel orientation angles for fed sheet 
network 

Panel edge lengths along spanwise cut for 

fed sheet network 

Geometry parameter X 

Geometry parameter V 

Total number of panel orientation angles 

for free sheet networks 

Total number of panel orientation angles 

for fed sheet networks 


XX 


F 

FX 

GG 

GX 

EX 


(1) - Potential 
(2,3,4) - Velocity vector 

(5) - Source strength 

(6) - Doublet strength 
(7,8,9) - Vorticity vector 
(10,11,12) - Unit upper surface normal n 
(13,14,15) - Control point z 
(16,17,18) - Upper surface normal n 
(19,20,21) - Normal vector to panel edge V 
(22,23,24) - Panel edge vector t 

Value of analysis boundary condition 
(ANLBC), or design boundary condition 
(DESBC), or edge matching condition (MATBC) 
Perturbation of boundary condition with 
respect to each of the variables XX 
Vector of panel force boundary condition 
(SFCBC and EGFBC) 

Perturbation of panel force boundary 
condition with respect to each of the 
variables XX 

(1,J) - Perturbation of force boundary 
condition fy with respect to Jth 
variable of XX 
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Common 

Block 


CM03 


CMOS 


CNTRL 


CNTRQ 


Variables 

Description 

(2,0) - Perturbation of force boundary 

condition f^ with respect to 0th 
variable XX 

NX 

24 

NX2 

48 

NX3 

72 

NTS IN 

Standard system file INPUT for card reading 

NTSOUT 

Standard system file OUTPUT for printing 

NSAV 

Disk file to save intermediated iteration 
results 

ERROR 

Error code from subroutine RHEAD or WHEAD 

I PASS 

Code for by passing the opening of random 
access file NSC2 

JPASS 

Code for by passing the opening of random 
access file NSC2 

ZC 

Control point "z 

ENC 

Upper surface normal ~n ^ 

ZNC 

Unit upper surface normal n 
n®Vf 

ZNCG 

ZCP 

Perturbation of control ^oint with respect 
to panel corner points P^ 

ENCP 

Perturbation o^r^ with respect to panel 
corner points P-j^ 

ZNCP 

Perturbation o;^ n with respect to panel 
corner points P^^ 

Perturbation of n®^f with respect to 

panel corner points Pj 

Normal vector to panel edge ^ 

ZNCGP 

AN 

EL 

Panel edge vectoj 9. 

ANP 

Perturbation of v with respect to panel 
corner points P __ 

ELP 

Perturbation of £ with respect to panel 
corner points P-\ 

ZSTC 

“+• 

^ara^te^ s and t of control point Qq, 
?s> ?b 3 ^t hyperboloidal 

panel of control point 

ICH 

Edge control point characterization 


= 0 Real 

= 1 Doublet value matching 
= 2 Doublet normal derivative matching 
= 3 Doublet tangential derivative matching 

I PC 

Panel no. of control point 

ISC 

Side no. of edge control point 

IZC 

Control point index along an edge 

JCN 

Control point index 
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Common 



Block 

Variables 

Description 


KC 

Network no. of control point 

COMPRS 

AMACH 

Mach number 


BETAMS 

1 - h2„ 


BETAM 

- m2co 


ALPC 

Angle of attack in degrees 


BETC 

Yaw angle in degrees 


COMPD 

Compressibility direction unit vector 


ACOMP 

Compressibil ity matrix for metric A 


BCOMP 

Compressibil ity matrix for metric B 


AROTC 

Orthogonal matrix transforming reference 
coordinates into wing axis coordinates 

CPLOT 

XMX 

Maximum X scale for the printer plot 


YMX 

Maximum Y scale for the printer plot 


XMN 

Minimum X scale for the printer plot 


YMN 

Minimum Y scale for the printer plot 


KFP 

Network no. of free sheet to be plotted 


NPLT 

Number for free sheet networks to be 
plotted _< 4 


ISY 

Symbols used in printer plot 

CRWI 

NCDQ 

Number of control point defining 
quantities per block 


NSC 

Number of control point defining quantity 
blocks in buffer 


NRC 

Current record in buffer 


NTC 

File on which control point defining 
quantity blocks are stored 

DMUDZ 

DMU 

Perturbation of panel singularity spline 
with respect to neighboring panel corner 
points 

DRW I 

NDDQ 

Number of singluairty spline derivatives 
per block 


NSD 

Number of singularity spline derivative 
blocks in buffer 


NRD 

Current record in buffer 


NTD 

File on which singularity spline 
derivative blocks are stored 

EDGIN 

KEDG 

Network of point which abuts given point 


MEDG 

Row index of point which abuts given point 


NEDG 

Column index of point which abuts given 
point 


lEGDA 

Cumulative index of points which abut 
other points 
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Common 

Block 

Variables 

Description 

EGMTCH 

ISDCHR 

Network edge control point characterization 


IRCHR 

=0 No control point 

= 1 to 4 Control point matches doublet 

strength along abutment to 
which side 1 to 4 belongs 
= 5 Control point forces doublet 

strength to vanish 
= -1 to -4 Control point matches 

vorticity parallel to edge 
along abutment to which side 
1 to 4 belongs 

Network corner control point 

FMCOF 

XREF 

characterization 
= 0 No control point 

= 1 to 4 Control point matches doublet 

strength along abutment to 
which side 1 to 4 belongs 
= 5 Control point forces doublet 

strength to vanish 
= -1 to -4 Control point matches 

vorticity parallel to edge 
along abutment to which side 
1 to 4 belongs 


YREF 

X, Y, Z coordinates of moment center 


ZREF 

SREF 

Configuration reference area 


BREF 

Reference span 


CREF 

Reference chords 


DREF 

Reference height 

INDEX 

NT 

Network type 


NM 

Network row number 


NN 

Network column number 


NZ 

Network mesh points number 


NP 

Network panels number 


NS 

Network singularity strength parameters 


NC 

number 

Network control points number 


NZA 

Accumulated network mesh points number 


NPA 

Accumulated network panel numbers 


NSA 

Accumulated network singularity strength 


NCA 

parameters number 

Accumulated network control points number 


NNETT 

Total number of networks 


NZMPT 

Total number of mesh points 
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Common 


Block 

Variables 

Description 


NPANT 

Total number of panels 


NSNGT 

Total number of singularity strength 
parameters 


NCTRT 

Total number of regular control points 


NCTRTE 

Total number of regular and special 
control points 

INDX 

INDX 

Key indices for random access file TAPE4 

LSINGC 

DSDFS 

Vector relating source strength o to 
neighboring singularity parameters 


DDDFS 

Matrix relating doublet strength )j and 
to neighboring singularity 

parameters 


DD 

MaJ;^rix relating doublet strength M and 

to the 9 canonical panel doublet 

parameters 

MSPNTS 

ZM 

Mesh points coordinates 

LSINGV 

SLV 

Vector consisting of source strength a , 
doublet strength y and n07y at a 
control point 


SLVP 

Perturbation of source strei^gth a , 
doublet strength y , and n®Vywith 
respect to panel corner points 

LSQSFC 

ZK 

X,Y,Z coordinates of mesh points used in 
least squares fit 


WTK 

Weights used in least squares fit 


AK 

Generalized inverse from least squares fit 


NO 

= 2 for quadratic fit (6 terms) 
< 2 for linear fit (3 terms) 


NPK 

Number of data points used in least 
squares fit 

NCONS 

PI 

7r 


PI2 

27T 


PHI 

1/47t 


AKAP 

47T 

NFAJ 

NF 

Number of singularity strength parameters 


NG 

Number of panel orientation angles 


NH 

Number of geometry parameters and 


NFG 

NF + NG 


NGH 

NG + NH 


NFGH 

NF + NG + NH 
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Common 



Block 

Variables 

Description 


NK 

Number of equations of twist function 


NGHK 

NG + NH + NK 


NKGHK 

NF + NG + NH + NK 

NITF 

JT 

Iteration No. 


ITMX 

Maximum number of iterations requested by 
the user 


JTP 

= 0 Printing of detail physical flow 
quantities 

= 1 Printing of type 1 or 2 force/moment 
data only 


ITPRIN 

Printing of detail physical flow 
quantities occurs at every ITPRIN iteration 


NDZA 

Type of iteration 

= 0 subiteration 

= 1 Iteration without force boundary 
condition 

= 2 Full iteration using Quasi-Newton 
scheme 

= 3 Full iteration using least squares 
method 


ITCODE 

Code for data check and type of iteration 
= 0 Data check 
= 1 Sub iteration 

= 2 Full iteration using Quasi-Newton 
scheme 

= 3 Full iteration using Least Squares 
method 


ITVRCP 

Printing code for variables, residuals and 
corrections of full iteration results 
= 1 Printout 
=0 No printout 

PANDQ 

CP 

Corner points of the given panel in 
reference cooj;;di nates 


PC 

^)» Qs’ Qst given panel 


AQ 

Transformation matrix from reference to 
panel near plane coordinates 


AQI 

Inverse of AQ 


DIAM 

Diameter of the given panel 


AST 

Matrix relating 9 canonical doublet 
or 3 linear source a to neighboring 
singularity parameters 


ITS 

Index array for neighboring singularity 
parameters 


INS 

Number of neighboring singularity 
parameters 
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Common 

Block 


PIVM 


PRNT 


Variables 

NCS 

ITS 

ICS 

IPN 

KP 

DVDS 

AMU 

DVDZ 

DVDP 


IGEOMP 

ISINGP 

ICONTP 

IBCONP 

ISINGS 

IPLOTP 

IPTIME 

IPRAIC 

IPAJAC 

IPNPIC 

IPSOLV 


E)escription 

Number of parameters (quadratic or linear 
coefficients) defining panel doublet or 
source distribution 
Panel singularity type 
= 1 Source 
= 2 Doublet 

Index of collasped side, and is equal to 0 

if no collapsed side 

Index of the given panel 

Network number of the given panel 

Potential and velocity vector at a control 
point induced by singularity distribution 
on a panel 

9 canonical doublet jj or 3 linear 
source o strength parameters 
Perturbation of potential and velocity 
vector with respect to the given control 
point 

Perturbation of potential and velocity 
vector with respect to the corner points 
of a panel 

All variables in this common block are 
printing codes for diagnostic purpose, 

= 1 Printout 
= 0 No printout 

Panel defining quantities 

Panel singularity distribution defining 

quantities 

Panel control points defining quantities 

Boundary conditions 

Singularity grid on each network 

Printer plot of cross sections of vortex 
sheet 

Elapsed CPU time for various programs and 
subrouti nes 

Influence coefficients for each control 
point 

Analytic Jacobian 

Far field and near field information 
Out-of-core equations solver information 
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Common 

Block 

Variables 

Description 

PRWI 

NPDQ 

NSP 

NRP 

NTP 

Number of panel defining quantities 
Number of panel defining quantity blocks 
in buffer 

Current record in buffer 

File on which panel defining quantity 

blocks are stored 

REAIC 

NRAIC 

= 0 All influence coefficients will be 
computed 

= 1 Part of the influence coefficients 
affected by the perturbation of 
geometry will be recomputed 

SKRCHS 

PANQ 

Buffers containing multiple blocks of 
panel and singularity distribution 
defining quantities 


CNTQ 

Buffers containing multiple blocks of 
control point defining quantities 


DMUQ 

Buffer containing temporary multiple 
blocks of panel defining quantities; 
buffer containing multiple blocks of 
singularity spline derivatives 

SKRCHl 

DUMSKl 

Scratch array 

SOLN 

S 

Singularity (doublet and/or source) 
strength 

SOLNT 

NAIC 

NRHS 

NANS 

NSCl 

NSC2 

NSC3 

NSCR 

Disk files (see File Usage) 


IRAY 

Array contaiing input specification for 


out-of-core equations solver 

(1) - Number of words in the scratch 


array 

(2) - I/O device for the coefficient 

matrix 

(3) - File argument for coefficent matri 

(4) - I/O device for the solution matrix 

(5) - File argument for solution matrix 

(6) - Scratch I/O device 

(7) - Scratch I/O device 


Common 



Block 

Variables 

Description 

(8) - Not used 

(9) - I/O matrix for the right-hand side 

matrix 

(10) - File argument for right-hand side 

matrix 


MTITLE 

Title of solution 

SRWI 

NSDQ 

Number of panel and singularity spline 
defining quantities 


NSS 

Number of panel and singularity spline 
defining quantity blocks in buffer 


NRS 

Current record in buffer 


NTS 

File on which panel and singularity spline 
defining quantity blocks are stored 

SYMM 

NSYMM 

= 1 for asymmetric about X-Z plane 
= 0 otherwise 

TFMQ 

FC 

Accumulated force coefficients of networks 


FMC 

Accumulated moment coefficients of networks 


TCA 

Accumulated surfaces area of networks 

UPDIND 

NUP 

Network update index 


NAT 

(1,K) - Sequence no. of network to which 
side 1 of Kth network is attached 
(2,K) - Sequence no. of network to which 
side 4 of Kth network is attached 


NSD 

(1,K) - Side no. of network to which side 
1 of Kth network is attached 
(2,K) - Side no. of network to which side 
4 of Kth network is attached 


NCR 

(1,K) - Leading corner point no. of 

network to which side 1 of Kth 
network is attached 
(2,K) - Leading corner point no. of 

network to which side 4 of Kth 
network is attached 

ZMD 

ZMTH 

Dependence of corner point coordinates on 
thetas 


ZMAL 

Dependence of corner point coordinates on 
lambdas and nus 


NZMTHA 

Cumulative number of thetas in each 
network 


NZMALM 

Cumulative number of lambdas and nus in 
each network 



7.5 Linkage Map of Overlay Programs and Subroutines 

In the following map, the subroutines followed with * are CDC system 
routines. 

(This linkage map is obtained from using a program CALLMAP written by 
Gary Bills of Boeing Computer Services Company.) 


7.5.1 MAP OF OVERLAY (MAIN. 0. 0) 


HAP OF 0VEPLAr(HAIN,0,0) 
ENTRY PT, A378V8 


E eof* 

OVERLAY* U.0),C2,0), C3,0) ,(6,0) 
OPENHS* 


fcWHEAD— 

^CTHNS — 

—STBNS — 

— ZERO 
— READMS* 


LOCATE' 


XRNS 


TRNS 


-C 


EOF* 

NIT* 

XFERA 

UNIT* 

XFERA 


p! 

P 


• CROSS 

• UVECTM- 

■CPNTCL--+'ZERO 

•CROSSM 
'CHAB 


■SORT* 


UlNlMU 


-SINFCC 

—XFERA 



SINCD 

CHAB 

SINCS 


l-ANLBC 


^CMAB 


l-DESBC 



ZERO 

XFERA 

CAMAB 

CHAB 



ZERO 

XFERA 

CAHAB 

CROSS 

CHAB 


E CAHAB 

MATBC*— -ZERO 
OVERLAY*, 15,0) 


A37BV8 


CLOSMS* 


VIP 



MAP 

ENTRY 


or OVERLAY(MAXN,Of 0) 
PT. A378VB 


f*UPDATE< 


•FUNC- 


XFERA 
COS* 
SIN* 
^OVERLAY* (3,0) 


►CTRNS* 


•FGCAL* 


■TRNS 


hCPNTCL 


^STRNS- 


"I 


CROSS 

UVECTK- 

ZERO 

CROSSM 

CMAB 


HSNGCAL— 


•TRNS 


• ZERO 


fc-SINFCC’ 


UNIT* 

XFERA 

■SORT* 

UNIT* 

XFERA 


•SINCD 
►CMAB 
► SXNCS 


HCMAB 
SINCD 


HXFERA 

HREADMS* 

kCHAB 


kANLBC* 


• DESBC* 


^-MATBC* 

[ ZERO 
CTRNS* 


►ZERO 

•XFERA 

►CAMAB 

•CMAB 

•ZERO 

•XFERA 

•CAMAB 

■CROSS 

►CMAB 

•ZERO 


•THNS« 


•UNIT* 


MAP or OVERLAY (MAIN,0,0) 
ENTRY PT« A378V8 


l-HCAL 


Hstrns 


TRNS 



XFERA 

UNIT* 

XFERA 


^CPNTCL 



CROSS 

UVECTM 

ZERO 

CROSSM 

CMAB 


—SORT* 


^ZERO 


hSINFCC 



SINCD 

CHAB 

SINCS 


Hsngcal 


-XFERA 

-READMS* 

-CMAB 


CMAB 

SINCD 


1-SFCBC 



DESBC 

CMAB 

ZERO 



ZERO 

XFERA 

CAMAB 

CROSS 

CMAB 


h-VADD 


^EPNTCL 



MNMOO 

UVECTM 


— SORTf 



ZERO 

XFERA 

CAMAB 

CROSS 

CROSSM 

CMAB 

SORT* 


^UNIT* 



HAP OF OVEHLAY(HAIN,0,0) 


ENTRY PT, A378V8 


#»STRNS< 


l-KCAL- 


*MAG* 
•SORT* 
•CMAB 


■TRNS* 


-I 


•CROSS 
•UVECTM- 

•CPNTCL-^-ZERO 

•CROSSM 
•CHAB 


•XFERA 


•SORT* 


•SORT* 


■VRCOUT VIP 

l-JTFLOM— 4*0VERLAY* (4/0), (5, 0| 

• WHEAD — LOCATE— EOF* 
•SORT* 

•XFERA 


f-UPDATE 



XFERA 

COS* 

SIN* 


UoVERLAY* (2,0) 


^CTRNS- 


|»CPNTCL 


HSTRNS- 


•TRNS 




CROSS 

UVECTH* 

ZERO 

CROSSH 

CMAB 


•TRNS 


h: 


UNIT* 

XFERA 

•SORT* 

UNIT* 

XFERA 


l-SNGCAL* 


ZERO 
h-SINFCC- 


t CHAB 
5INC0 


•SINCD 

•CMAB 

•SYNCS 


t-XFERA 



MAP 

ENTRY 


OF OVERLAY(MAIN»0,0) 
PT. A378V8 


^FUNC- 


■FGCAL- 


Hhcal- 


•READHS* 

•CHAB 


|-ANLBC< 


Hoesbc 

■HATBC- 

•ZERO 

•CTRNS- 

l-STRNS' 


-i 


•ZERO 

•XFERA 

■CAMAB 

•CMAB 

ZERO 

XFERA 

CAMAB 

CROSS 

CMAB 

ZERO 


■TRNS 


>TRNS 


"I 


•CROSS 
'UVECTH' 

•CPNTCL— ^ZERO 

►CROSSM 
•CMAB 


UNIT* 

XFERA 

UNIT* 

XFERA 

•SORT* 


f-ZERO 


hSINFCC- 


I-sngcal-h 


•SINCD 

•CMAB 

•SINCS 


VCMAB 
SINCD 


•XFERA 

•READMS* 

•CMAB 


'-ii 


• ZERO 
•XFERA 
•DESBC— +-CAMAB 
•CROSS 



HkP OF OVERLAY(MAIN,0,0) 
ENTRY PT. A378V8 


^-KCAL• 


hSFCBC 


hVAOD 


kEPNTCL 


l-EGFBC 


^HAG 


•STRNS' 


— C 


4-CMAB 


GRAB 

ZERO 


|-HNMOD 
l—UVECTM— — 


SORT* 



l-CPNTCti 


•MAG« 

t SQRT* 
CMAB 


CROSS 

CROSSM 

CMAB 

SORT* 


■TRNS 


-c 


-i 


CROSS 

UVECTM- 

ZERO 

CROSSN 

CMAB 

'SORT* 


UNIT* 

XFERA 

'SORT* 


t-lTFLOW- 


-VIP 

-VRCOUT— — VIP 

-OVERLAY* (4,0), (5,0), (6,0) 


hWHEAD- 


hRHEAD* 


—LOCATE' 

C LOCATE* 
EOF* 


'EOF* 

■EOF* 


t XFERA 
SORT* 


^UPDATE 


-Ei 


XFERA 
COS* 
SIN* 
OVERLAY*(2,0) 



MAP 

ENTRY 


OF OVERLAY(MA1N,0,0) 
PT. A376V8 


^FUNC> 


■FGCAL 


#-CTRNS 


TRNS 


UNIT* 

XFERA 


h-CPNTCL 



CROSS 

UVECTH 

ZERO 

CHOSSM 

CHA0 


—SORT* 


f-STRNS— TRNS 



UNIT* 

XFERA 


^ZERO 


l-SINFCC 



SINCD 

CMAB 

SINGS 


l-SNGCAL 



CHAB 

SINCD 


E XFERA 

REAOMS* 

CMAB 



ZERO 

XFERA 

CAMAB 

CHAB 

ZERO 

XFERA 

CAHAB 

CROSS 

CHAB 

ZERO 


f»ZERO 


t-CTRNS— — TRNS 



UNIT* 

XFERA 


t-STRNS 


TRNS 



UNIT* 



MAP 

ENTRlf 


OF OVCRLAY(MAIN,0,0) 
PT. A378V8 


t-HCAL- 


kCPNTCL 


-i 


CROSS 

UVECTH* 

ZERO 

CROS5M 

CHAB 


4«XFERA 


•SORT* 


^ZERO 


l-SINFCC 


hSNGCAL-H 


-E 


SINCD 

CMAB 

SINCS 


kxFERA 

REAOMS* 

hCMAB 


j-CMAB 

*-SINCD 


hSFCBC 


h-VADD 


4 DESB( 

CHAB 

ZERO 


P 

DESBC — 1» 


ZERO 

XFERA 

CAMAB 

CROSS 

CHAB 


|— MNHOD 

-EPNTCL I 

^-UVECTH^^ - 


SORT* 


h-EGFBC- 


Hhag* 


l-STRNS- 


^ZERO 

-XFRRA 

-CAMAB 

—CROSS 
— CR05SM 

4 — CMAB 
SORT* 


►TRNS 


-C 


UNIT* 

XFERA 


ACROSS 


MAP OF OVERLAY(HAXN,0,0) 
ENTRY PT, A378VB 


1-LSFliOW* 


KVIP 

|-VRCOUT— — VIP 

J-OVERLAY* (4,0), (5,0), (6,0) 


-Lctrs -E 


Hvips 


LkcAIi- 


^CPNTCL 


i 


MAG< 
SORT* 
V-CHAB 


UVECTH« 

ZERO 

CROSSM 

GRAB 

•SORT* 


•SORT* 


hRHEAD 


kWHEAD- 


-E 


h«TDCOM 


XFERA 

h-HTBSUB' 


-C 


LOCATE— 

EOF* 

LOCATE— 

ZERO 

WRITMS* 

PEADKS* 

VIP 

SORT* 


•EOF* 


■EOF* 


•VIP 
■ VIP 


*-XFERA 

♦“UPDATE— 4-COS* 

I l-SIN* 

l“OVERLAY*f 7,0) 


r-CTRNS« 


•TRNS 


— I 


_£ 

-CPNTCL— ^ 


CROSS 

UVECTM- 

ZERO 

CROSSM 

CMAB 


UNIT* 

XFERA 

•SORT* 


-►-UNIt* 



MAP 

ENTRY 


OF OVERLAY(MAIN,0»0) 
PT. A378V8 


|*STRNS< 


•TRNS' 


>XFERA 


l-FUNC- 


ff-ZERO 


hSlNFCC- 


•FGCAL' 


l-SNGCAL— I 


l-XFERA 

■READHS* 

•CMAB 


t-ANLBC 

•DESBC 

•MATBC- 
•ZERO 
•CTRNS* 

HSTRNS' 


•SINCD 

•CMAB 

■SINGS 


t CMAB 
SINCD 


-E 

-i 


ZERO 

XFERA 

CAMAB 

GMAB 

ZERO 

XFERA 

CAHAB 

CROSS 

CMAB 

ZERO 


•TRNS 


•TRNS 


-c 


■-I 


•CROSS 
■UVECTM' 

•CPNTCL— +-ZERO 

•CROSSM 
'CMAB 


UNIT* 

XFERA 

UNIT* 

XFERA 

■SORT* 


|-ZERO 


•►■SINCD 


133 



HAP OF UVERI.AY(KAIN,0,0) 
ENTRY PT, A378V8 
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7.5.2 MAP OF OVERLAY (DATA. 1. 0) 


MAP OF OVERLAlf(DATA,l,0) 

entry pt, input 


•SORT* 

'OUAONT 

'GADNET 


< 


ZERO 

SORT* 


I ipXFERA 

^VORTEX I SHECEN-— »SQRT* 
I ^TRWAKE--- XFERA 


HTRWAKC' 


Kcamber« 


•XFERA 
:CEI 
fAKE^ 

•XFERA 


p-ISCAL- 
XFERA 


h-INTPl 


^IHTP2- 


•PIDENT> 


'BIMSCH 


•BINSCH 


k-CONCAM 



ASIN* 

COS* 

SORT* 


INPUT 


^COS* 

ksjN* 


hROTATE 



COS* 

SIN* 


^SCALE2— — ALOGIO* 
UVSPLOT' ---SUBPON 


►-•btcal 



SORT* 


SORT* 



f 


MAP or OVCRLAIT(OATA,1,0] 
ENTRY PT. INPUT 
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7.5.3 MAP OF OVERLAY (NETGCSj 1 . 0) 


HAP UF OVERLAY(NETGCS,2,0): 

ENTRY PT. CONFIG 

CONFIG— -OVERLAY* ( 2 r 1 > r ( 2, 2) , ( 2,3) 
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7,5.ti MAP OF OVERUW (NETGCS. 


2 . 


1 ) 


MAP OP OVERLAY(NETGCSr2»l)l 
ENTRY PT, TGEOHC 


♦-GEOMC- 


•PCCAL 

■VADD 

■COMPIP 

■SORT* 


t-SURFIT 


t-ISCAL- 




^IPTRNS* 


•CROSS 


•MAG" ' ■ 

—SORT* 

SORT* 


JORDAN 


PIDENT-- 

-^QRT* 


^XFERA 

ITRNS — 

1 


^•UNIT* 


TGEOMC- 


^•CSTPRT' 


•SECOND* 



7.5.5 MAP OF OVERLAY (NETGCS^ 1 . 2) 

MAP OF UVEHLAY(hFTGCS,2»2): 
entry PT, TSiNG 


•GCPCAL 

'GRDIND — PI DENT— SORT* 


Hptrns- 


C XFERA 


•UNIT* 


•SING* 


•tPROJ' 


•UNiPANi 

■VADO 


•CHAB 

•UVECTm^m 

■SORT* 

'PROJ 

-VADD 

-CRAB 

H C n D X A 


VADD 


•CMAB 



p-ZERO 


■LSOSFX 


h-tSOSF- 


•SCMPKT* 


•ISTRNS' 


|-*PSINTP— 

l-V.P 

UoCBHiX -c 
L.UTPS 


SORT* 

t-pi dent* 


r XFERA 
VADO 


OCBHIX 
VIPS 

•SORT* 


•VIP 

SORT* 

VIP 


•ITRNS 


-c 


UNIT* 

XFCRA 


-ZERO 

-GCPCAL 

-MNSWCH 

-GROIND— PI DENT— SORT* 
— XFERA 
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MAP 

OF 

OVERLAY(NETGCS,2,2) 1 

ENTRY 

PT. 

TSING 


i-PTRNS 


TRNS 



UNIT* 


^DASPL 


-MNMOD 

-INCHPR 

-XFEHA 


r-ZERO 


koiSTNC 



MAG- 

VAUD 


SORT* 


t P I DENT— SORT* 
EDGLS 


l-BLCCAL 


MSROTM 



TRANS 

UVECT 

CROSS 


SORT* 


TRANS 

MNHOD 


P CHAB 
UVECT 


h-LPROJ 


PROJ 



t MAG- 
VADD 


r UNlPAN«**CHAB 
rZERO 


SORT* 

VADO 

CRAB 

SORT* 



PSINTP 

VIP 

DCBHTX 

VIPS 


—VIP 



SORT* 

VIP 


Vlsosf 


t SORT* 

PIDENT*— SORT* 



0 VERLAY(NETGCS, 2 , 2 ) 


HKP 

ENTRY 


OF 

PT, TSINC 


l-DDSPL* 


hBLCAL- 


j“ 

"C 

t-lSTRNS^— 'URNS -C 


kSCHPKT- 


•ZEHO 

>MNMOD 

'XFERA 

•VADD 


•GRDIND* 


kPTKNS- 


hMNMOO 

XFERA 


■PIDENT' 


'TRNS 


-C 


UNIT* 

XFERA 

SORT* 

XFERA 

UNIT* 


^ZERO 


•LSQSFX 


LSOSF— 

[ UNIPAN' 

LPROJ- 

'MNMOD 


t-BLDCAL' 


pPSINTP 

J-V.P 

LdCBHTX— 

WviPS 


t SQRT* 
PIDENT- 
— CMAB 

|-CHAB 
>-UVECT — 


►PROJ 


MAG* 
l-VADD 


DCBHTX 
VIPS 

•SORT* 

■SORT* 

VAUD 

CMAB 

SORT* 


•VIP 

SORT* 

VIP 


MAP 

ENTRY 


or OVERLAY(NETGCS,3»2)t 
PT, TSINC 


h-SCMPKT 


h-TRANS 


h-MSROTM 



TRANS 

UVECX 

CROSS 


SQRT« 


r-ZEKO 


h-Lsosrx 


J-LSOSF— 

tuNIPAN' 

-LPROJ — 


HEHDSPL' 


f-PSINTP— 

l-VIP 

j-DCBHTX— 

Lvips 


kSORT* 
i-PIDENT- 
CHAB 

^CMAB 

VUVECT — 


•PROJ 


t MAG — 
VADO 


l-HNMOD 
TRANS 


hMSROTH 


Hdcdgls 


-DISTNC 


-E 

**C 


TRANS 

UVECT- 

CROSS 


tPIDENT- 


VADD 

SORT* 


•ZERO 

• XFERA 

• VADD 


DCBHTX 
VIPS 

•SORT* 

•SORT* 

VADD 

CHAB 

SOR"^ 


•VIP 

SORT* 

VIP 


•SORT* 


•SORT* 



MAP OF 0VERLA¥(NETGCS,2,2) : 
ENTRY PT. TSING 


TSING-' 


L -T' 

i-ISTRNS — ITRNS— ^ 


•GCPCAL 

>GRDINO< 


•PIDENT— SORT* 


>DSDSPL< 


*BL5CAL< 


LSQSFX* 


■UNIPAN — CMA8 


CMAB 

UVECT- 


■PSINTP-— VIP 
•VIP 


•DCBHTX' 


SORT* 

PIDENT-— SORT* 


'-C 


pTRANS 
•MSROTM— l-UVECT- 
i-CROSS 


143 



NAP 

ENTRY 


OF OVERLAY(NETGCSf 2»2)l 
PT. T51NG 


r-DEDGLS 


HCNESPL' 


I-schpkt-H 


'ZERO 

•XFERA 

•VADD 


UlSQSF' 


Hunipan< 


hLPROJ 


hMNHOD 

TRANS 


h-MSROTN- 


^XFERA 

Hedgls 


•ZERO 


h-LSOSFX 


i 


PSINTP' 

VIP 


t SORT* 
PIDENT- 

*CMAB 



CMAB 

UVECT' 


PROJ 


— C 


MAG— 

VADD 


•TRANS 

• UVECT* 

• CROSS 


DCBHTX 

VIPS 

•SORT* 

• SORT* 
VADD 
CMAB 
SORT* 




■VIP 

SORT* 

VIP 


•SORT* 


hDISTNC 


t-MAG* 

UvADI 


■SORT* 


kPIDENT* 


VADD 

SORT* 


•-UNIT* 


MAP 

ENTRY 


or OVERLAY(NETGCS,2r2): 
PT. TSING 


t-DwSPL* 


>1STRNS — — ITRNS 


l-EDGLS 

Vdei^gls 


•SCMPKT 


< 


XFERA 


VADD 


XPERA 


XFERA 

PIDENT 

——SORT* 


pXFERA 

PTRNS- 

— trns— 

i-UNIT* 

HNHOD 

ZERO 


D1* A 

DISTNC 

pMAG 

LvADD 



»-lSTRNS — --1TRNS 


-c 


IWIT# 

XFERA 


WCSTPRT- ■■SECOND* 



7.5.6 MAP OF OVERLAY (NETGCS. 2. 3) 


MAP OF OVERLAYCNeTGCS»2,3)S 


ENTRY PT. TCNTRL 


TCNTRL- 


rCONTRL-H 


•GCPCAL 

•GRDIND—— PI DENT— -SORT* 


■ STRNS — 

■XFERA 

■ICTRNS- 


•TRNS 


T 


•ITRNS 




•CROSS 
•UVECTM* 

•CPNTCL-^-ZERO 

•CROSSM 
•CHAB 


UNIT* 

XFERA 

XFERA 

UNIT* 

•SORT* 


t STRNS* 
XFERA 


•CONFRC-H 


•TRNS 


— c 


klCTRNS^— ITRNS 


UNIT* 

XFERA 

XFERA 

UNIT* 


"l! 


CROSS 

UVECTM— SORT* 
UCPNTCL— f-ZERO 

CROSSM 
CHAB 


WCSTPRI— SECOND* 
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7.5.7 MAP OF OVERLAY (AICGEM, 3. 0) 

MAP GIF 0VEWLAYtAICGEN,3,0)S 
ENTRY PT. VINFCC 


r-ZERO 


l-FFHPMG 



DQGM 

ASITJO 

RRZAB 

ZERO 


t CTRNS — 

XFERA 
READHS* 

h-STRNS-— 


VINFCC—^ 


>TRNS* 


■TRNS« 


♦-EIFLU' 


•XFERA 

•UNIT* 

■XFERA 

■UNIT* 

■CPTLS 

•VADD 

Hxxadj 


•HNMOD 

■CMAB 


l-ZHPROJ« 


-i 

-c 


MAC* 


■SORT* 


VADD 


MAG* 


■SORT* 


•PIDENT* 
.KAG— — 


VADD 

•SORT* 
•SQR 7* 


■DISTNC 


Kxfera 


J-MAG* 

^VADD 


•SORT* 


CMAB 

M-CPOSSM 

-cpntcl— i-zero 

-UVECTH— — SORT* 
l-CROSS 


l-MNMOD 


'SORT* 

•VADD 

■VMUb 

■CMAB 



MAP 

ENTRY 


or OVERLAY(AlCGEN«3rO)t 
PI, VINFCC 


h>NRPTED 




SORT* 

xfera 

COMPIP 

VADD 


HINSIPE 



CMAB 

CROSS 

XFERA 

VADO 


t UVECT — SORT* 
CROSS 


kEIVC 


hCPNTCL 



CHAB 

CROSSH 

ZERO 

UVECTH 

CROSS 


— "SORT* 


-SINFCC 

-EPNTCL 


“C 


SINCS 

CHAB 

SINCD 

UVECTH 

HNHOO 


■-ZERO 


—SORT* 


I— CPNTCL 



CHAB 

CROSSH 

ZERO 

UVECTH 

CROSS 


kZERO 

Ktypcal— HNMOD 


SORT* 


l-PIVC 


INFLU 



CSCAL- 

CMAB 

SORT* 

COMPIP 


VADD 


SCROSS 

ACROSS 

AGAT 

CAXAB 

RRAABT 



MAP 

OF 

0VEKLAY(AICGEN»3,0) > 

ENTRY 

PT. 

VINFCC 


^CAMAB 


^FFPIC 


-XFKKA 

-RRZ.ABT 

-CXAB 

-RRZATB 

-COHPIP 

-VAOD 

-SORT* 

-CHAB 

-CROSS 

-CSCALl 

-ZERO 


hPlFCAL 


ONEPIF 


r-SUBSBX 



IPArj-^- 

LqcOF— — — 

l-CMAB 


ALOG* 

ATAN2* 

SORT* 

ZERO 

CSA8 

ZERO 


Kbl-anq 



VADD 

UVECT 

CMAB 

SORT* 

TRANS 


E UVECT— — SORT* 
CROSS 
PCCAL 


L -CAMAB 
-CROSSH 
-CROSS 

tilNVOR— -ECAL 
—A LOG* 
-SORT* 
-CHAB 


t CAMAB 

VRITH5* 


SORT* 



HAP 

ENTRY 


OF OVERLAYCAICGEN, 9»0}l 
PT, VINFCC 


|-AICFOR< 


•CTRNS' 

^XFERA 

i-ZERO 

• STBNS’ 


HPIVC* 


■TRNS 


*TRNS 


‘CPNTCL 


XFERA 

UNIT* 

XFERA 

UNIT* 


p-CMAB 

UCROSSM 
— 4-ZERO 

^UVECTH«— 

Across 


SORT* 


•ZERO 

•TYPCAL' 


•MNMOD 


^CSCAL- 
^CHAB 

l-SQRT* 

^COHPI 


•INFLU 


•VADD 


•FFPIC- 


•SCROSS 

•ACROSS 

•AGAT 

•CAXAB 

•RRAABT 

•CARAB 

•XFERA 

■RRZAHT 

•CXAB 

■RRZAT8 

•CORPIP 

•VADD 

•SORT* 

■CMAB 

•CROSS 

•CSCALl 

•ZERO 


ALOG* 

ATAN2* 


UpIFCAL' 


•ONEPIF 




MAP 

OF 

OVERLAY (AlCGEN, 3,0)! 

ENTRY 

PT. 

VINFCC 


hPLANO 



VADD 

UVECT 

CMAB 

SORT* 

TRANS 


SORT* 


E UVECT— — SORT* 
CROSS 
PCCAL 


l-LINVOR 



CAMAB 

CROSSM 

CROSS 

CCAL 

ALOG* 

SORT* 

CMAB 


t CAMAB 
WRITMS* 

-CSTPRT— SECOND* 
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7.5.8 MAP OF OVERLAY (JACGEN. A. 0) 

MAP OF OVERLAyCjACGCN, 4,0)1 
ENTRY PT. AJGEN 


l«Mi 


>ZERO 
_ -KRZAB 

•FFHPKG* 

" .ASITJQ 
>DUGM 
.XHEAD—— LOCATE' 


k*ZTHET< 


■EOF* 


^COS* 

VsiN* 


■ OAST< 


-GRDIND— PIDENT — 
-PIDENT— SORT* 

—SORT* 

pUNlT* 

-STRNS' TRNS — 

Lxfera 

-CMAB 

-MNMOD 


-XFERA 

HZERO 


^VADD 
-DISTNC— 1 

Umag— 

—SORT* 


^ZERQ 


^PIDENT- 


h-EMDSPL- 


■DISTNC 


kDEDGLS 


■HSROTM 


TRANS 

h-HNMOD 


“C 

-E 


SORT* 

VADD 

RAG— —SORT* 


CROSS 

UVECT' 

TRANS 


4 VADD 
HAG — 

PROJ» 


■SORT* 




SORT* 

CHAB 

VADD 
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MAP 

OF 

0VERLAY(JACGEN,4,0) 

ENTRY 

PT. 

AJGEN 


fc 


h-UNlPAN- 


*-tSOSF- 


UVECT- 

CWAB 

•CHAB 


r-PIDEWT- 

hSQRT* 


-SORT* 


• SORT* 


rvips 


r-DCBHTX 


•L505FX* 


VZERO 


-c 


VIP 


SORT* 


■VIP 

•PSINTP— — VIP 


r-BLDCAL-H 


■M5ROTK 

-TRANS 

-MNNOD 


-E 


CROSS 

UVECT- 

TRANS 


K-LPHOJ- 


l-UNIPAN- 


■-LSOSF- 


p-VRDO 
HKAG — 


•SORT* 


-proj 


-c 


UVECT~ 
^CMAB 

•CHAB 

r-PIOENT' 

HSQRT* 


i-LSOSFX 


i 


SORT* 

CMAB 

VAOO 

-SORT* 

•SORT* 

VIPS 

DCBHTX 




VIP 

SORT* 
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HAP 

ENTRY 


OP OVERLAY(JACGEN«4,0>f 
PT. AJGEN 


t ZERO 
IDTRNS 


l-CTRNS' 


>TRNS 


•DASTPS-H 


IDTRNS 


ZERO 

XPERA 

PIDENT— 


DISTNC 

OEDGLS 


“C 


ITRNS 


-c 


SORT* 
VADD 
MAG — 

XFERA 

UNIT* 


Lzero 


SORT* 


ITRNS 



XPERA 

UNIT* 


UNIT* 

XPERA 


C VIP 

PSINTP — VIP 


rZERO 


hCPNTCL 



CROSS 

UVECTM 

ZERO 

CROSSM 

CMA8 


—SORT* 


Udvjcdt 


t STRNS— TRNS 
CHAB 


UNIT* 

XPERA 



CROSS 

UVECTM 

ZERO 

CKOSSH 

CMAB 


—SORT* 
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MAP 

OF 

OVERLAY (JACGEN, 4 rO) 1 

ENTRY 

PT. 

AJGEN 


t-EIFLU 


t-XF£RA 


l-DISTNC 




t MAG — 
PIOENT- 


LxXADJ— m 


Hvadd 


Lcptls -c 


VADD 

MAG- 

SORT* 

SORT* 

VADD 

HAG— 

VADD 

HAG- 


SORT* 


SORT* 


SORT* 


-HNHOD 

-CHAB 

F CROSS 

UVECT— SORT* 


l-ZHPROJ 


h-lNSIDE 


h-WRPTED 



VADD 

XFERA 

CROSS 

CHAB 

VADD 

CDMPIP 

XFERA 

SORT* 


-CHAB 

-VHUL 

-VADD 

-SORT* 


I-DEIVC 


l-CPNTCL 



CROSS 

UVECTH 

ZERO 

CROSSH 

CHAB 

ZERO 


—SORT* 


HAP 

OP 

UVERLAYCJACGEN,4^0)k 

ENTRY 

PT, 

AJGEN 


-SNGCAL' 

-EPNTCL 

^ZERO 

^DISTNC 

UzERO 

UtYPCAL’ 


•SINFCC 


■CMAB 

•SINCD 


-E 


SINCD 

CHAB 

5INCS 


pHNHDD 

Luvectk— - 


SORT* 


•VADD 

•HAG- 


SORT* 


•INFLU 


-E 


• MNHOD 

COMPIP 

SORT* 

CMAB 
CSCAL — 


•VADD 

^PCCAL 

Across 

UUVECT- 


l-PLANQ 


kPIFCAL- 


•ONEPIF-H 


-CMAB 

-QCOF— 

l-UNIPAN- 


■SORT* 


f-OBLlQU 

^CSCAL- 

tcHAB 
L-sort* 


-f 


TRANS 

SORT* 

CHAB 

UVECT« 

VAOD 


•SUHSBI 


•ZERO 
•CMAB 

•—ZERO 

^SQHT* 

UaTAN2* 

kAi.nr.* 


SORT* 
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MAP 

ENTRY 


OF OVERLAY(JACGEN,4,0)t 
PT. AJGEN 


hDPlVC 


-CMAB 

-XFERA 

-PCCAL 


I-FFPIC 


-ZERO 

-C.SCALl 

-CROSS 

-CMAB 

-SORT# 

-VADO 

-CQMPIP 

-RPZATB 

-CXAB 

-RRZABT 

-XFERA 

— CAMAB 

-RRAABT 

-CAXAB 

-AGAT 

-ACROSS 

-SCROSS 


fc-LINVOR 



CHAB 

SORT* 

ALOG* 

ECAL 

CROSS 

CROSSM 

CAMAB 


^DTRNS 


TRNS 


UNIT* 

XFERA 


■-DFDTH 



MNHOO 

DZDTH 

CAMAB 



HNMOD 

XFERA 

ZERO 


4-UNlI* 


MAP or OVERLAY (JACGCN,4^0)I 
ENTRY PT. AJGEN 


fSTRNS' ■ " TRNS 


UdTRNS'* ™ 'TRNS 



XFERA 

UNIT* 

XFERA 


l-CPNTCL 



CROSS 

UVECTM 

ZERO 

CROSSH 

CMAB 


—SORT* 


AJGEN 


UXFERA 

hzERO 


l-SNGCAL 


t READMS* 
CMAB 


ZERO 


SINFCC 



SINCD 

CMAB 

S1NC5 


CMAB 

SINCD 



ZERO 

XFERA 

CAMAB 

CMAB 

CMAB 

CAMAB 


DFDTH 



MNMOO 

DZOTH 

CAMAB 


MNHOD 

XFERA 

ZERO 


MAP 

entry 


OF OVERLAY(JACOEN,4»0)I 
PT. AJGEN 


l-DESBC 


h-MATBC- 


-I 


ZERO 

XFERA 

CAHAB 

CROSS 

CMAB 


ZERO 

[ ZERO 
CTRNS- 


■TRNS 


-c 


UNIT* 

XFERA 


h-DVJCOT- 


r-ZERO 


hCPNTCL 

-STRNS- 

•CMAB 


■| 


CROSS 

UVECTM- 

ZERO 

CROS5H 

CMAB 


■TRNS 


-c 


•SORT* 

UNIT* 

XFERA 


i-EIFLU- 


rMNMOD 


hCPNTCL 


h-XFERA 


hOISTNC 


-i 


CROSS 

UVECTM- 

ZERO 

CROSSM 

CMAB 


■SQRT4 


HI 


VADD 

hag— 


'SORT* 


•MAG* 


HPIDENT’ 


LxXADj_j^ 


•SORT* 

•SORT* 

VADD 

MAG— 


•SORT* 


VADD 



HAP 

ENTRY 


OF QVERLAYCJACGENt 4«0) t 
PT. AJGEN 



■CPTLS 


VADD 
HAG — 


r CROSS 

UVeCT ■ — SORT* 


k-INSIOE 


^-NRPTED 



VADO 

XFERA 

CROSS 

CMAB 

VADD 

COMPIP 

XFERA 

SORT* 


^CHAB 

l-VHUL 

^VADD 

L-SORT* 



CROSS 

UVECTM 

ZERO 

CRUSSH 

CHAH 


SORT* 


r»ZEHO 


hSINFCC 



SINCD 

CMAB 

SINCS 


t CMAR 
SINCD 


-c 


HNMOD 

IJVECTH« 


•SORT* 


■^VADD 


SORT* 
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II 


MAP 

OF 

OVERLAY (JACGEN^4«0) t 

ENTRY 

PT. 

AJGEN 


-DISTNC 

“ZERO 

-TYPCAL 



SORT* 


— MNMOD 


hiNFLU 



COMPIP 

SORT* 


CKAB 

CSCAL 


VADD 


E PCCAL 
CROSS 

UVECT— SORT* 


l-PLANQ 


|“0BLIQU 

l-CSCAb- 

I-cmab 

l-SQRT* 


CMAB 

SORT* 



TRANS 

SORT* 

CMAB 

UVECT- 

VADD 


•SORT* 


I-PIFCAL 


ONEPIF 


CMAB 

OCOF— — ZERO 
UNIPAN— — CMAB 


►-SUBSB1 



ZERO 

SORT* 

ATAN2* 

AtOG* 


kDPIVC 


“CMAB 

“XFERA 

“PCCAL 


■“FFPiq 


“ZERO 

“CSCAL 1 

“CROSS 

“CMAB 

“SORT* 

“VADD 

“COMPIP 

“RRZATB 

“CXAB 
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I 


MAP 


ENTRY 


OF OVERLAY(JACGEN«4fO)l 
PT. AJGEN 


•RRZABT 

•XFERA 

>CAMAB 

>RRAABT 

•CAXAB 

-AGAT 

•ACROSS 

•SCROSS 


CHAB 


kDTRNS" 


1-LlNVOR 


'TRNS 



^DFOTH 



NNMOO 

XFERA 

ZERO 


hSTRNS TRNS 


L.DTRNS— TRNS 



UNIT* 

XFERA 

UNIT* 

XFERA 


l-CPNTCL 



CROSS 

UVECTH 

ZERO 

CROSSM 

CMAB 


—SORT* 


♦-ZERO 


162 


HAP 

ENTRY 


OF 0VERLAY(JACGEN, 4#0)t 
PT. AJGEN 


l-DFRC 


h-SNGC* 


L 



SINFCC 


CMAH 

SINCD 



SINCD 

CMAB 

SINGS 


h-XFERA 


I^REAONS* 

HCHAB 


t*SFCBC 



DESBC 


CHAB 

ZERO 



ZERO 

XFERA 

CAHAB 

CROSS 

CHAB 


h»VADD 


h-FDER 



CHAB 

CAHAB 

DFDTH 



HNHOD 

DZOTH 

CAHAB 



HNHOD 

XFERA 

ZERO 


l-ePNtCL 



HNHOD 

UVECTM 


—SORT* 


K-egfbc 



ZERO 

XFERA 

CAHAB 

CROSS 

CROSSM 

CHAB 



HAP 

ENTR» 


OF 0VERLAY(JACGEN,4,0)| 
PT. AJGEN 


L-MAG 


SORT* 


fSTRNS 


TRNS 



UNIT* 

XFERA 


l-CPNTCL 



CROSS 

UVECTM 

ZERO 

CROSSM 

CHAB 


— SORT* 


HOKCAL 


.mag — 

-SORT* 

-CHAB 

-ZERO 

-CAMAB 


SORT* 


l-DEDTH 



MNMOD 

DZDTK 

CAMAB 


MNMOD 

XFERA 

ZERO 


^CSTPRT — SECOND* 



7.5.9 MAP OF OVERLAY (SOLVER. 5. 0) 

MAP OF OVFRLAY(SOLVFR,S,0): 
entry PT. TINVER 


TINVER 


pRHEAD -c 

•OPENMS# 

■ZERO 


EOF* 

LOCATE- "EOF* 


■BLOCKR 


r-DlRBLK 


< 


-BLKFCT 

— openms* 

-ZERO 


RVEC— 

WRITMS* 


EOF* 


hCHOPR 


^MURGR 


-c 

-c 


RVEC-’ — "EOF* 
WRITMS* 

READMS* 

WRITMS* 


P FILSET 
REPSOL 


h-PTDCDM 


pREADMB' 

kcSHAB 

Uaramax 

PRbSUB' 

h-wRITMB' 


READMS* 


VIPS 

WRITMS* 


hPCDCOH 


l-PFSUB 


~L 

-c 


XCHR 

SLIN 

XCHR 

VIPS 


l-PTSOLV 


h-PTFSUti 


P READMB— RE ADMS* 
CSMAB 


I-PF5DB 



^WRITMB 


XCHR 

VIPS 

WRITMS* 


t README— READMS* 
CSMAB 
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HAP OF OVERLAY (SOLVER,5tO) 
ENTRY PT, TINVER 


t RBLOCK 
CSTPFT 


KTBSUB 



BSUB ■ 
WRITMB— — ■ 


VIPS 

WRITHS* 


— PRPIV — ■ ■-READHB— READMS* 

»WHEAD ■■ I ■ ■'LOCATE' -EOF* 
■READHB- ■ -READMS* 

-WVEC 


SECOND* 



7.5.10 MAP OF OVERLAY (RESULT.' 6. 0) 


MAP OF aVERL.AY(ReSUl.T,6,0)t 
ENTRY PT, OUTPUT 


-'HNHOD 

-STRNS- 

-XFERA 

-CPNTCL 


t UVECT- 
CROSS 


h-SNGCAL-H 


OUTPUT* 


'TRNS 


-I 


CROSS 

UVECTM. 

ZERO 

CROSSM 

CMAB 

-SORT* 


UNIT* 

XFERA 

•SORT* 


•ZERO 


h-SINFCC 


-£ 


SINCD 

CMAB 

SINCS 


CMAB 
WSINCD 


l-CMAB 

SORT* 

j-ze«o 

-CTRNS — 

— READMS* 
^CSCAL — 

-CPCAL -E 


•TRNS 


-c 


•VADD 


SORT* 

CSCAL' 

CMAB 


UNIT* 

XFERA 

•VADD 


ZERO 

^STRNS- 

^CROSS 

VUVECT- 

Lubliou 


•TRNS 


— c 


UNIT* 

XFERA 


-E 


•SORT* 

TRANS 

SORT* 

CMAB 

UVECT. 


•SORT* 



HAP 

ENTRY 


OF DVERLAY(RESULT,6,0) t 
PT, OUTPUT 


|-FMCAL< 


•UNIPAN 
•CCAL 


CKAB 

ZERO 

ECAL 


-CSCAL VADD 

-CHAB 

-SORT* 


kLPROJ 


>VADD 

■MAG' 


■PROJ 


•UVECT' 

•CHAB 


SORT* 

CMAB 

VADD 

SORT* 


Hlsosf 


P PIDENT 
SORT* 


Hlsqsfx 


t t ^ZERO 

l-VADD 

VSPLOT— SU8P0N 
CSTPRT— SECOND* 


—SORT* 


I— VIPS 


l-DCBHTX 



VIP 

SORT* 


t VlP 

PSINTP— -VIP 



REFERENCES 


1. Mask ell, E. C., “Some Recent Developments in the Study of Edge Vortices," 
Proceedings of 3rd Congress of Int. Counc. Aero, Sci., 1962, pp 737-749, 
Spartan Books, Inc., Washington, 1964. 

2. Peckham, D. H., "Low-Speed Wind-Tunnel Tests on a Series of Uncambered 
Slender Pointed Wings with Sharp Edges," RM 3186, British Aeronautical 
Research Council, 1961. 

3. Smith, J. H. B,, "Improved Calculations of Leading-Edge Separation from 
Slender Delta Wings," RAE Tech, report 66070, March 1966. 

4. Tinoco, E. N. and Yoshihara, H., "Subcritical Drag Minimization of Highly 
Swept Wings with Leading Edge Vortices," Paper No. 26, A6ARD-CP-247 , 
January 1979. 

5. Wentz, W. H., "Effects of Leading Edge Camber on Low-Speed Characteristics 
of Slender Delta Wings," NASA CR-2002, 1972. 

6. Kuhlman, J. M., "Analytical Studies of Separated Vortex Flow on Highly 
Swept Wings," NASA CR-3022, November 1978. 


169 





2. Governmeni AcceKion No 


3 Recipient's CdialOQ No 


1 Report No 

NASA CR-3279 


15. Supplementary Notes 

Langley Technical Monitors: James M. Luckring and Ward E. Schoonover, Jr. 

Topical Report 

16 Abstract 

An improved panel method for the solution of three dimensional flow about wing 
and wing-body combinations with leading edge vortex separation is presented. 
The method employs a three-dimensional inviscid flow model in which the con- 
figuration, the rolled-up vortex sheets, and the wake are represented by 
quadratic doublet distributions. The strength of the singularity distribu- 
tion as well as shape and position of the vortex spirals are computed in an 
iterative fashion starting with an assumed initial sheet geometry. The meth- 
od calculates forces and moments as well as detail surface pressure distribu- 
tions. Improvements include the implementation of improved panel numerics 
for the purpose of eliminating the highly non-l-inear effects of ring vortices 
around doublet panel edges, and the development of a least squares procedure 
for damping vortex sheet geometry update instabilities. 

The documentation is divided up into two parts: 

Volume I - Theory Document 

Volume II - User's Guide and Programmer's Guide 

Volume II contains instructions for the proper set up and input of a problem 
into the computer code. Program input formats and output are described. A 
description of the computer program and its overlay structure is also present- 
ed. 

18, Distribution Statement 

Unclassified-Unl i mi ted 

Subject Category 02 

f this page) 21, No ot Pages 

id 173 

•for sate by the National Technical Information Service. Springfield. Virginia 22161 


22 Puce' 
$ 8.00 


1 7 Key Words (Suggested by Author(s) I 

Leading Edge Vortex 
Panel Method 
Three-Dimensional 
Separation 

20. Security Classif. 

Unclassif 


19. Security Classif. (of this report) 

Unclassified 


4. T itle and Subtitle 

An Improved Panel Method for the Solution of Three- 
Dimensional Leading- Edge Vortex Flows. Volume II - 
User's Guide and Programmer's Document 

7. Author(s) 

E. N. Tinoco, P. Lu, and F. T. Johnson 


9 Performing Organization Name and Address 

Boeing Aerospace Company 

P. 0. Box 3999 

Seattle, Washington 98124 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


5 Report Date 

July 1980 


6 Performing Oigam/ation Code 


8, Performing Organization Report No 


10 Work Unit No 


11 Contract or Grant No 


NASI -151 69, NASI -15 


11 T ype of Report and Period Covered 

Contractor Report 
December 


14. Sponsoring Agency Code 


Jk 


NASA-Langley, 1980 











